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OPENING  REMARKS 


Wayne  Howell 

Chief  of  Fire  Safety  Branch 
FAA  Technical  Center 

Good  morning,  my  name  is  Wayne  Howell  and  I  am  Chief  of  the  Fire 
Safety  Branch  here  at  the  Federal  Aviation  Administration  Technical 
Center.  I  would  like  to  welcome  you  to  the  first  aircraft  fire  math¬ 
ematical  model  workshop.  As  most  of  you  know,  the  FAA  has  been  in¬ 
volved  in  aircraft  fire  math  modeling  for  a  few  years  and  has  achieved 
some  accomplishments.  We  thought  it  would  be  good  to  invite  some  ex¬ 
perts  in  fire  modeling,  but  not  necessarily  just  those  involved  in  the 
aircraft  fire  modeling  area,  to  have  a  technical  exchange  of  informa¬ 
tion  today  and  tomorrow.  We  would  like  to  show  you  what  we  have 
accomplished  and  what  we  are  doing.  Hopefully  you  will  learn  some¬ 
thing  on  aircraft  fire  safety  RfiD  work  and  possibly  we  will  also 
improve  our  program  as  a  result  of  some  of  your  critique  and  comments. 
It  is  a  very  informal  conference  and  we  would  like  you  to  relax  and 
enjoy  the  presentations.  When  you  make  a  comment,  I  would  like  to  ask 
you  to  please  stand  up,  speak  up  a  little  bit  louder,  and  identify 
yourself  so  that  your  comments  can  be  recorded. 

I  would  like  to  explain  to  you  how  the  aircraft  fire  math  modeling 
work  relates  to  the  overall  activities  at  the  FAA  Technical  Center. 

Some  of  you  are  familiar  with  the  Technical  Center's  operations  and 
some  are  not,  so  I  would  like  to  start  off  by  showing  you  first  of  all 
that  the  FAA  Technical  Center  (Figure  1)  is  the  most  extensive  proving 
ground  of  aviation  safety  systems  in  the  United  States.  The  Technical 
Center's  Mission  is  shown  in  Figure  2.  It  also  has  international 
recognition  because  many  FAA  regulations  formed  from  technical  data 
developed  here  are  used  as  international  standards.  Particularly,  we 
are  the  leaders  in  the  field  of  aviation  safety  standards.  The  re¬ 
search,  development,  and  testing  that  we  do  here  at  the  Center  evolves 
into  new  concepts,  new  procedures  in  communications,  navigation,  air 
traffic  control,  and  aircraft  and  airport  safety. 
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The  Technical  Center  is  located  about  12  miles  from  Atlantic  City 
and  has  approximately  1500  employees  (Figure  3) .  The  Center  goes 
back  to  about  1958.  Prior  to  1958,  there  was  a  technical  center  in 
Indianapolis  under  the  old  Civil  Aeronautics  Administration.  When 
the  FAA  Act  was  established  in  1958,  this  Center  was  set  up  here  in 
New  Jersey  in  place  of  a  former  naval  air  station.  The  Center  has 
5,000  acres  and  about  1,000,000  square  feet  of  building  space.  The 
new  Technical  and  Administration  Building  has  500,000  square  feet  of 
floor  space  and  houses  close  to  1,000  people. 

The  Center  has  the  most  modern  airport  in  the  United  States  with 
a  10,500  foot  long  runway.  The  newest  most  advanced  aviation  concepts 
are  being  tested  here. 

In  order  to  test  out  new  concepts  in  communication/navigation/ 
aircraft  safety  and  air  traffic  control,  the  Center  has  a  complete 
cross  section  of  aircraft  from  a  helicopter  and  small  propeller  type 
airplane  up  to  a  large  jet  Boeing  727  shown  in  Figure  4,  We  are  very 
well  equipped  here  to  perform  our  mission. 

The  Technical  Center  organization  chart  is  shown  in  Figure  5. 

Mr.  Joseph  Del  Balzo  is  the  Director.  There  are  four  divisions  which 
do  the  actual  research/development  and  test/evaluation  work.  The 
Systems  Test  and  Evaluation  Division  test  and  evaluate  air  traffic 
navigation  and  communications  procedures  and  facilities.  The  Systems 
Simulation  and  Analysis  Division  simulates  air  traffic  control  pattern 
or  configurations.  They  can  simulate  any  air  traffic  control  pattern 
or  configuration  in  the  world.  They  have  simulated  air  traffic  patterns 
of  Chicago  O'Hare  Airport,  one  of  the  largest  and  busiest  airports  in 
the  world.  New  procedures  and  new  techniques  for  more  efficiently 
handling  the  air  traffic  at  large  airports  are  being  developed.  The 
Aircraft  Safety  Development  Division,  which  I  will  go  into  in  more 
detail  later,  is  where  the  fire  modeling  work  is  being  done.  The 
Airport  Technology  Division  is  looking  at  approach  and  runway  config¬ 
urations.  The  Center  also  has  some  tenant  organizations.  The  Flight 
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ORGANIZATION 


Figure 


Inspection  Division  here  at  the  Center  has  six  or  seven  jet  aircraft 
that  flight  inspect  the  communications  and  navigation  facilities  in 
the  eastern  part  of  the  United  States  to  determine  if  they  are 
functioning  properly  and  accurately. 

Before  the  systems  are  implemented,  they  are  tested  and  evaluated 
here  at  the  Center,  particularly  in  the  air  traffic  control  and 
conununicat ions/navigation  areas.  In  the  aircraft  safety  area,  our 
chief  product  is  technical  criteria  which  are  the  basis  for  new  regu¬ 
lations  or  revisions  of  current  regulations. 

The  Aircraft  Safety  Development  Division  is  the  only  division 
here  at  the  Center  which  has  the  complete  overall  responsibilities  for 
research,  development,  test,  and  evaluation.  At  this  point,  I  would 
like  to  introduce  Dr.  Roy  Reichenbach,  the  Division  Chief.  The  Air¬ 
craft  Safety  RfiD  Program  is  shown  in  Figure  6.  If  you  have  any 
questions  concerning  this  division's  operations,  Dr.  Reichenbach  is 
certainly  available  to  answer  those  questions.  The  division  consists 
of  a  Propulsion  and  Fuel  Safety  Branch,  Crashworthiness  Branch,  Opera¬ 
tions  Branch,  and  Five  Safety  Branch.  The  Fuel  Safety  Branch  has  kaD 
work  going  on  in  antimisting  fuel,  which  is  designed  to  reduce  the 
post-crash  fire  hazard.  Antimisting  fuel  is  a  fuel  which  has  been 
modified  by  adding  a  polymer,  In  a  crash  situation,  the  fuel  spills 
out  of  the  fuel  tank  and  atomizes  to  flammable,  small  droplets.  The 
polymer  added  to  the  fuel  prevents  it  from  becoming  small  droplets 
and  thereby  reduces  the  fire  hazard.  Another  way  of  trying  to  prevent 
a  fire,  of  course,  is  to  design  the  airplane  to  withstand  a  certain 
crash  impact.  The  Crash  Worthiness  Branch's  program  is  to  develop 
and  strengthen  the  fuselage  and  fuel  tanks  to  withstand  higher  impacts. 
Under  flight  safety,  the  Operation's  Branch  is  looking  at  the  airplane 
itself,  trying  to  design  the  airplane  to  be  more  compatible  with  the 
pilot  to  reduce  pilot  error. 

The  Aircraft  Safety  Development  Division  has  approximately 
$15,000,000  worth  of  facilities  at  the  Center  (Figure  7).  These 


8 


£AA  AIRCRAFT  SAFETY  R&D 


<  o 
co  2 
< 

cu  £ 

rv  00 


W  * 
®  H  CO 
* ^NH 
ik  h  m  i 

_  w  •  Q 

2b.HU 
M  <  <0  | 
B  W  h  i 

U  £ 


25 

5  5 

as  c5 

LU  XT 

*  oo 


<5 

CO  d 


-<  LU 

cr  S 

O  Q. 

as  o 

M  I 

<c  LU 

£ 

o 


£5 

85 

5  as 

£  03 

CO 

< 

or 

o 


o  (0 

gl  88 

siisv 

wp<<N 
J  Cl  06  <0 
to.  X  3  m 

to) 

0.  O 


tfj 

<A  v0 
,  U  I 
f-  2  O  ® 

to.  M  M  rl 

<  x  in  i 
epic 

O  X  CM  u 
X  O  ®  i 
~  < 

<  s  < 

<  ■** 


O 
Q  2 
2  < 
<  a: 
oo 


2fc 

CO  i-u 
—I 

ZD  < 
n  lO 

O 

as  _i 

Q_  CU 


Ik  O  B 

CM  .H 

a  m  i 

u  i  a 


O  Z  <30 
NNPIh 

v)  x  m  i 
Jl*  Ifl 
p  X  CM  Cl] 
C  C  S  I 
-2-< 
?  X  < 


FIGURE  7.  AIRCRAFT  SAFETY  DEVELOPMENT  DIVISION  TEST  FACII TTIES 


facilities  were  designed  to  simulate  environmental  conditions  like  post¬ 
crash  fire,  in-flight  fire,  a  crash  situation  or  many  of  the  kind  of 
hazards  we  are  trying  to  reduce  or  prevent.  We  have  a  five-foot  wind 
tunnel  in  which  we  can  run  tests  on  small  jet  engines  to  determine  the 
fire  detection  and  extinguishing  requirements  for  the  engine.  We  have 
an  engine  test  facility  in  which  we  are  able  to  provide  air  into  the 
cells  to  simulate  flow  through  the  engine  cell  itself  which  simulates 
airflow  conditions  while  the  airplane  is  in  flight.  We  have  a  drop 
test  rig  in  which  we  can  determine  ways  and  means  of  containing  fuel 
in  the  wing.  We  have  a  catapult  and  track  facility  and  several  fire 
test  facilities.  We  have  a  component  laboratory  in  which  we  do  lab¬ 
oratory  fire  tests. 

The  Fire  Safety  Branch's  mission  is  illustrated  in  Figure  8.  The 
major  mission  is  to  improve  and  develop  fire  safety  standards  for  air¬ 
craft.  In  addition,  the  branch  is  developing  fire  protection  systems 
for  the  Air  Force. 

The  aircraft  fire  safety  work  covers  propulsion  systems,  fuel 
systems,  airframes,  cabin  related  components,  and  airport  fire  fight¬ 
ing  systems.  I  would  like  to  point  out  the  fact  that  we  actually 
design  our  own  unique  fire  test  facilities.  The  engineers  establish 
the  specifications  and  work  with  the  architect  to  insure  that  the  unique 
requirements  for  fire  testing  are  included. 

A  typical  example  is  the  full-scale  fire  test  facility  which  was 
just  completed  (Figure  9).  This  is  the  largest  in-door,  full-scale 
fire  facility  operated  by  the  federal  government.  It  is  185  feet  long, 
75  feet  wide  and  45  feet  high.  It  is  capable  of  housing  a  wide  body 
jet  inside,  with  the  wings  and  the  upper  tail  cut  off.  Currently,  we 
use  a  surplus  C-133  fuselage  and  an  8'xlO'  pool  fire  outside  the  fuse¬ 
lage  to  simulate  a  wide  body  jet  postcrash  fire  situation. 

A  new  chemistry  laboratory  (Figure  10)  is  under  construction  which 
will  be  completed  by  September  1981  and  it  will  be  utilized  to  study 
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toxicity  aspects  of  the  materials  in  a  postcrash  fire  situation. 

In  the  middle  of  the  airport,  there  is  a  200  foot  diameter  burn  pit 
in  which  we  are  able  to  dump  10,000  gallons  of  jet  fuel,  set  fire  to 
it  and  test  the  effectiveness  of  different  extinguishing  systems 
and  foam  agents.  We  have  tested  both  dry  chemical  and  wet  foam  to 
extinguish  an  aircraft  fire.  A  dual  agent  rapid  response  system  is 
very  effective  in  putting  out  the  fire  in  a  very  short  time  period 
and  is  affordable  by  small  aviation  airports. 

The  Technical  Center  has  been  recognized  for  its  expertise  and 
facility  capabilities  for  approximately  35  years  in  the  aircraft  fire 
safety  work.  As  a  result  of  this  expertise  and  facilities,  the  Air¬ 
craft  Fire  Safety  Branch  has  an  active  fire  protection  program  for 
the  Air  Force.  We  have  sclved  a  lot  of  fire  protection  problems  in 
military  aircraft.  The  major  efforts  are  outlined  in  Figure  11. 

The  Air  Force  has  lost  about  six  F-lll  aircrafts  due  to  in-flight 
fires  and  has  requested  the  FAA  to  investigate  the  cause.  An  F-lll 
fuselage  is  placed  on  a  fire  test  pad  without  the  wings  and  tail  sur¬ 
faces.  The  airflow  supplied  from  a  nearby  jet  engine  compressor 
bleed  air  into  the  engine  intake  and  simulates  actual  airflow  con¬ 
ditions  that  the  F-lll  would  have  during  in-flight  conditions.  A 
fire  is  set  within  the  engine  bay  to  test  various  new  fire  detection 
and  extinguishing  techniques  to  improve  the  fire  worthiness  of  this 
particular  aircraft. 

The  major  thrust  of  the  Aircraft  Fire  Safety  Branch  is  the  cabin 
fire  safety  progr~m  shown  in  Figure  12.  It  is  a  very  comprehensive 
program.  There  are  five  different  major  areas.  One  major  area  is 
Survival  and  Evacuation,  which  includes  management  of  people  during 
evacuation  and  applying  survival  aids  to  assist  people  in  getting 
out  of  the  aircraft.  The  second  major  area  is  laboratory  test  meth¬ 
odology  development.  The  objective  is  to  use  laboratory  test  methods 
for  material  selection.  The  material  selected  should  be  less  flammable, 
produce  less  smoke  and  be  less  toxic.  The  third  major  area  is  fire 


16 


3- 


1 VI 


management  and  suppression  which  is  to  develop  effective  techniques 
of  managing  firs.  The  fourth  major  area  is  standards  and  improvements 
which  includes  risk  analyses,  data  bank  and  specific  material  improve¬ 
ments.  The  fifth  area  is  postcrash  cabin  fire  hazard  characterization 
where  fire  math  modeling  is  being  accomplished.  It  is  important  to 
characterize  the  pcstcrash  fire.  The  main  characteristics  measured 
are  fire  progression,  gas  temperatures,  radiation,  and  convective  heat 
fluxes.  Full-scale  fire  tests  and  math  modeling  are  conducted  in 
parallel  to  achieve  maximum  results. 

The  final  recommendations  from  our  fire  program  are  expected  to 
be  formalized  in  August  1983.  The  standards,  improvements,  and  accept¬ 
ability  criteria  will  be  ready  by  that  time.  The  aircraft  fire  safety 
program  is  moving  fairly  well  and  we  hope  this  symposium  today  will 
accelerate  it.  That  briefly  gives  you  the  scope  and  the  relevance  of 
the  fire  math  modeling  program. 

I  would  like  to  make  a  few  administrative  announcements.  Tomor¬ 
row  morning  we  will  have  a  tour  of  our  facilities.  We  realize  that 
some  of  you  have  already  seer,  our  facilities  and  probably  will  not  be 
interested  in  joining  the  tour.  For  those  who  do  not  go  on  the  tour, 

I  have  reserved  two  small  conference  rooms  up  on  the  fourth  floor  in 
which  wc  can  continue  our  discussions. 

The  Math  Model  Advisory  Panel  will  meet  in  the  tower  room  at 
Resorts  International,  Wednesday  evening  at  7:30-10:00  p.m.  The 
minutes  of  this  symposium  will  be  summarized  and  there  will  be  a 
proceedings  published. 
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CONSTANTINE  3ARK0S 

Program  Manager,  Cabin  Fire  Safety  Program, 
Fire  Safety  Branch,  FAA  Technical  Center; 
w . j .  Mechanical  Engineering,  Rutgers.  Gus 
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AIRCRAFT  FIRE  SCENARIOS 


Constantine  (Gus)  Sarkos 
Fire  Safety  Branch 
FAA  Technical  Center 

Thank  you  very  much,  Wayne.  Good  morning.  The  aircraft  safety 
record  has  been  relatively  good  over  the  years,  especially  compared 
to  other  modes  of  transportation.  There  have  been  major  fluctuations 
in  the  statistical  data.  For  example,  last  year  there  was  only  one 
fatal  accident  which  resulted  in  fourteen  fatalities,  and  the  record 
was  excellent.  Yet  several  years  ago,  over  500  died  at  the  Canary 
Islands  during  the  collision  of  two  747s,  and  many  of  these  fatalities 
were  due  to  fire.  The  potential  exists  for  fire  in  a  commercial  trans¬ 
port  due,  of  course,  to  the  large  quantities  of  jet  fuel  carried  on 
board  and  also  due  to  the  cabin  interior  materials  which  are  of  interest 
to  the  workshop  on  fire  modeling. 

The  cabin  interior  of  a  modern  jet  transport,  in  this  case  a  wide- 
body,  is  obviously  far  different  from  that  of  a  residential  room  which 
has  been  the  type  of  enclosure  that  has  been  most  often  modeled  in  the 
past.  Typically,  the  carpet  is  of  wool  pile  construction.  Seats  are 
constructed  of  a  treated  urethane  foam  and  wool  nylon  upholstery  fabric. 
Polycarbonate  is  used  on  window  dust  screens  and  on  passenger  service 
units.  The  sidewall  panels,  ceiling  and  partitions  are  of  composite 
panel  type  construction.  The  length/width  ratio  of  an  airplane  cabin 
is  on  the  order  of  10  versus  one  to  three  for  a  residential  enclosure. 
The  materials  found  within  the  home  are  relatively  untreated  whereas 
those  used  in  airplanes  are  very  thermally  resistant  and  are  compliant 
with  the  FAA  flammability  regulation  which  is  basically  a  Bunsen  burner 
type  of  test.  Moreover,  the  materials  are  tested  by  the  airplane 
manufacturers  for  low  smoke  emissions  and  also  for  f ireworthiness  on 
the  basis  of  both  small  scale  and  large  scale  tests.  Many  of  these 
panels  have  flame  spread  indices  of  less  than  ten. 
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The  multilayer  panels  constitute  most  of  the  surface  area  within 
the  cabin  interior,  especially  in  the  upper  levels  of  the  cabin.  From 
a  fire  safety  viewpoint,  panels  are  clearly  the  most  important  materials 
system.  They  are  also  very  complex  in  terms  of  geometry  and  composition. 
A  typical  multiple  layer  panel  construction  is  shown  in  Figure  1.  It 
has  10  components  of  various  thermal  properties  within  the  same  panel 
construction.  A  composite  panel  used  on  the  three  types  of  wide  body 
jets  typically  consists  of  a  polyvinyl  fluoride  decorative  film  finish, 
fiber  glass  facings  impregnated  with  either  epoxy  or  phenolic  resin 
and  a  Nomex  paper  honeycomb  core.  The  behavior  of  a  composite  panel 
under  fire  almost  defies  description,  especially  compared  with  the  type 
of  materials  that  have  been  most  often  studied  in  the  modeling  work  by 
many  fire  researchers.  For  example,  plexiglass  is  slow  burning  and 
well  behaved  compared  to  composite  panels.  In  the  older  airplanes, 
vinyl  on  aluminum  was  used  on  sidewall  panels  and  vinyl  on  cloth  was 
used  for  the  upper  cabin  areas. 

There  are  three  types  of  aircraft  cabin  fire  scenarios,  as  shown 
in  Figure  2,  They  are  ramp  fires  (fire  occurs  usually  when  the  airplane 
is  unattended),  in-flight  fires,  and  postcrash  fires. 

A  ramp  fire  is  a  smoldering,  long  duration  type  of  fire  which  is 
more  like  a  residential  fire.  This  problem  is  not  addressed  in  the 
FAA  modeling  work. 

The  in-flight  fire  is  also  not  of  immediate  interest  to  the  current 
FAA  mathematical  fire  modeling  program.  For  in-flight  fire  safety,  we 
are  more  concerned  with  developing  an  early  fire  detection  system,  and 
an  effective  extinguishing  system  to  contain  the  fire. 

The  FAA  mathematical  fire  modeling  program  focuses  on  the  post¬ 
crash  cabin  fire,  which  accounts  for  all  the  fatalities  attributable 
to  fire  in  U.S.  air  carrier  accidents. 

The  aircraft  interior  arrangement  is  designed  for  rapid  evacuation. 
There  is  an  FAA  regulatory  requirement  that  complete  evacuation  be 
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accomplished  within  a  time  period  of  90  seconds  or  less  with  half  of 
the  exit  doors  open.  This  is  a  requirement  which  must  be  demonstrated 
by  the  airframe  manufacturers  and  airlines  during  a  live  evacuation 
drill.  Here  exists  another  major  difference  between  the  residential 
and  aircraft  fire  modeling  work.  The  time  of  interest  in  aircraft 
fires  is  0  to  5  minutes  while  in  a  residential  fire,  it  is  in  the 
neighborhood  of  an  hour  or  even  longer, 

I  would  like  to  show  a  sequence  of  slides  which  are  perhaps  the 
most  detailed  ever  taken  of  an  aircraft  accident.  The  accident 
occurred  several  years  age,  following  tire  blow-out  and  landing  gear 
collapse.  There  was  penetration  of  the  fuel  tank  and  a  major  fire 
erupted.  The  fire  was  on  the  left-hand  side  of  the  airplane.  The 
wind  which  was  blowing  right  to  left  had  a  very  important  bearing  on 
the  development  of  the  fire.  In  this  instance,  the  R-2  slide  (evacua¬ 
tion  slide)  caught  fire  while  the  R-l  slide  was  still  being  used.  The 
R-l  slide  failed  later  due  to  radiant  heat.  This  started  us  on  a  pro¬ 
gram  to  improve  the  heat  resistance  of  slides.  Now  the  airframe  man¬ 
ufacturers  are  devoting  attention  to  improving  the  heat  resistance  cf 
aircraft  evacuation  slide  materials.  In  this  case,  the  firefighters 
became  aware  of  the  accident  even  before  the  airplane  came  to  a  hair 

and  were  at  the  accident  site  within  100  seconds  after  the  initial 

ignition.  It  is  surprising  how  severe  this  fire  was  yet  there  were 
only  two  fatalities  probably  caused  by  disorientation  of  two  elderly 
passengers.  The  aircraft  fire  was  extinguished  in  an  estimated  time 
period  of  about  three  minutes.  The  fire  was  attacked  from  the  right- 

hand  side  while  the  larger  fire  was  actually  on  the  left.  The  orienta¬ 

tion  of  the  airplane  and  the  final  landing  resulted  in  an  open  space 
beneath  the  airplane.  The  fire  on  the  other  side  can  be  seen  through 
the  opening.  From  examination  of  past  accidents,  we  can  come  up  with 
three  major  characteristics  of  a  survivable  postcrash  fire.  The 
first  is  a  large  external  fuel  fire.  Practically  all  crash  accidents 
with  fire  involve  spillage  of  jet  fuel  although  there  are  a  few 
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exceptions.  If  there  is  a  major  breakage  of  the  fuselage  or  even  a 
separation,  the  burning  fuel  rather  than  the  involvement  of  interior 
materials  becomes  the  predominant  hazard.  Consequently,  in  a  realis¬ 
tic  fire  test,  there  should  be  an  opening  placed  in  the  vicinity  of 
the  fire  to  allow  flames  and  heat  to  penetrate  and  ignite  the  interior 
materials.  In  the  experiments  performed  by  the  FAA,  a  typical  door 
opening  adjacent  to  a  test  fire  is  used. 

The  next  question  is  how  to  treat  this  large  fire  adjacent  to  a 
long  airplane  fuselage.  A  lot  of  test  work  has  been  done  on  this 
subject  area.  One  treatment  to  this  problem  is  to  study  the  fire 
penetration  through  one  large  opening  and  basically  ignore  any  pene¬ 
tration  through  the  remaining  part  of  the  fuselage.  This  treatment 
is  valid  for  a  short  time  interval  in  a  wide-body  jet  which  is  con¬ 
structed  with  highly  fire  resistant  materials.  The  excellent  fire 
resistance  to  burn-through  penetration  of  the  interior  panel  con¬ 
struction  used  in  a  wide-body  jet  was  demonstrated  in  a  fire  acci¬ 
dent.  After  three  min"*es  or  longer  of  exposure  to  a  major  fire, 
there  was  significant  melting  of  the  aluminum  skin  but  no  flame 
penetration  to  the  interior. 

Pool  fires  have  been  studied  extensively  over  the  years,  The 

2 

radiant  heat  flux  is  relatively  invariant  at  about  14  Btu  ft  /sec 
for  pool  fires  of  three  feet  in  diameter  or  greater. 

2 

The  convective  flux  is  much  smaller  at  1  to  3  Btu  ft  /sec,  and 
is  dependent  on  the  size  of  the  fire.  A  plot  of  the  radiant  flux 
by  a  fire  plume  is  shown  in  Figure  3,  Assuming  the  fire  could  be 
treated  as  a  black  body  radiant  sphere,  the  receiving  heat  flux  at 
various  distances  are  calculated.  An  inverse  square  relationship 
for  the  decrease  in  radiation  versus  distance  is  obtained.  The 
practical  deduction  here  is  that  in  order  to  have  any  smoldering  or 
flaming  combustion  on  the  cabin  interior,  the  fire  has  to  be  adjacent 
to  the  fuselage.  A  fire  adjacent  to  an  opening  will  produce  very 
intense  radiant  heat  and  ignite  the  cabin  inside  materials.  The  flame 
penetration  through  the  opening  depends  on  wind  speed  and  direction 
and  location  of  other  openings. 
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The  FAA  has  studied  the  penetration  of  fire  through  fuselage 
openings  using  subscale  models  as  well  as  surplus  airplanes.  A  four- 
foot  diameter  fuselage  model  was  made  from  an  open-ended  cylinder  and 
an  opening  was  placed  on  the  side.  A  fuel  pan  was  placed  adjacent  to 
the  opening  to  simulate  the  postcrash  fire  scenario. 

The  FA\  also  studied  fire  penetration  using  a  surplus  DC-7  adja¬ 
cent  to  a  2C-?oot  square  fuel  fire.  There  was  a  fire  penetrating 
through  the  opening.  The  amount  of  the  penetration  is  dependent  upon 
the  wind  velocity  vector  as  well  as  the  placement  of  openings  away 
from  the  fire. 

The  ceiling  temperatures  inside  the  DC-7  fuselage  versus  time 
for  a  number  of  high  wind  cases  are  shown  in  Figure  4;  fire  was  up¬ 
stream  of  the  fuselage.  The  worst  condition  (high  temperature  in  a 
short  time)  occurred  when  the  downwind  door  was  open  and  the  upstream 
door  was  closed.  This  apparently  was  due  to  the  low  pressure  area 
created  by  the  wind  flow  over  the  aircraft  cylinder  creating  a  draft 
which  induced  flame  penetration  into  the  cabin.  Contrast  this  with 
a  case  where  the  upwind  door  was  open  and  the  downwind  door  was 
closed.  A  very  moderate  penetration  of  flame  and  resultant  low  build¬ 
up  of  heat  inside  the  interior  was  measured. 

Wind  in  the  aircraft  postcrash  fire  can  be  a  detrimental  factor 
to  hazard  development,  Wind  induced  flame  penetration  will  also  in¬ 
crease  radiant  heat  flux.  This  is  illustrated  in  Figure  5,  Radiant 
heat  flux  on  the  symmetry  plane  against  time  was  measured  for  the 
calm  wind  and  mild  wind  cases.  A  reasonable  agreement  on  heat  fluxes 
was  achieved  between  the  modeling  and  calm  wind  results  with  all  doors 
closed.  For  a  fluctuating  wind,  shown  by  the  dashed  curve,  the  radiant 
heat  fluctuated  above  the  calm  wind  pattern.  As  a  result  of  flame 
penetration,  the  radiant  and  convective  heat  fluxes,  smoke  and  gases 
inside  the  cabin  increased  as  the  fire  penetrated  further  into  the 
interior. 
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FIGURE  4.  DC 7  CEILING  TEMPERATURE  HISTORIES 


4  10  MPH  WIND 


FIGURE  5.  DC 7  SYMMETRY  PLANE  HEAT  FLUX 


The  bulk  of  the  FAA  full-scale  test  work  is  now  conducted  with  a 
C-133  wide  body  test  vehicle.  An  8'xlO’  fuel  pan  which  produces 
approximately  70  or  80  percent  of  the  radiant  heat  from  a  very  large 
pool  fire  is  placed  adjacent  to  a  door  opening.  The  initial  tests  in 
the  C-133  were  baseline  tests  to  determine  the  hazards  from  the  fuel 
fire  alone.  The  cabin  interior  was  bare.  We  found  out  that  from  an 
external  fuel  fire,  heat  and  smoke  were  much  more  hazardous  than 
carbon  monoxide.  We  did  not  measure  much  carbon  monoxide  inside  the 
airplane  from  a  fuel  fire. 

The  impact  of  the  external  fuel  fire  on  the  cabin  interior 
materials  is  of  greatest  concern  to  the  FAA  fire  safety  program. 

Seat  materials  are  chosen  as  initial  candidates  for  the  study.  We 
are  studying  seat  fires  resulting  from  ignition  by  a  large  pool  fire 
penetrating  through  an  opening.  We  are  looking  at  improved  cushion 
materials  which  will  be  a  viable  replacement  for  the  currently  treated 
urethane.  However,  about  a  year  or  more  ago,  we  did  run  one  test  in 
the  C-133  with  a  20-foot  section  furnished  and  lined  with  seats  and 
materials  (used  in  a  wide-body  jet)  that  were  provided  by  airframe 
manufacturers  and  various  suppliers.  Coincidentally,  we  are  running 
a  test  similar  to  this  today.  It  is  our  first  test  of  this  nature 
inside  our  new  fire  test  facility.  The  seat  cushions  are  protected 
by  a  fire  blocking  layer. 

A  test  was  set  up  to  illustrate  that  a  major  fuel  fire,  external 
to  the  airplane,  would  ignite  internal  materials  that  in  turn  would 
affect  the  passengers  survivability.  Basically,  information  pertain¬ 
ing  to  the  development  of  fire,  the  mechanisms  of  fire  development 
inside  the  aircraft  cabin,  and  the  buildup  of  hazards  were  collected. 

The  results  were  quite  obvious.  There  was  extensive  damage  near 
the  fire  door.  The  fuel  fire  did  ignite  the  interior  materials. 

There  was  fire  development  which  preceded  very  gradually  in  the  be¬ 
ginning  but  then  became  much  more  intense.  The  estimated  time  for 
survivability  in  the  cabin  in  this  particular  test  was  about  three 
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minutes.  There  would  have  been  virtually  no  hazards  from  the  fuel 
fire  alone.  The  hazards  were  strictly  due  to  the  involvement  of  the 
interior  materials. 

During  the  fire  test,  the  seat  next  to  the  door  ignited  very 
early  and  burned  rapidly.  However,  there  was  very  little  ignition  of 
other  materials  in  the  airplane  at  that  time.  One  minute  later  the 
seats  immediately  forward  and  aft  of  the  seat  in  the  opening  ignited. 
These  were  the  only  materials  which  were  burning  for  most  of  the  test. 
The  heat  produced  by  the  seat  was  rising  and  hot  gas  was  building  up 
at  the  ceiling.  This  caused  the  ceiling  panels  to  pyrolyze  and  dis¬ 
tort.  When  these  distorted  ceiling  panels  collapsed  onto  the  remaining 
seats  away  from  the  door,  a  very  rapid  growth  in  the  fire  was  observed. 
Tn  Figure  6,  temperature  measurements  versus  time  at  26  feet  aft  of 
the  fire  on  the  symmetry  plane  for  a  series  of  thermocouples  one-foot 
apart  from  the  floor  to  the  ceiling  were  plotted. 

It  is  very  interesting  and  perhaps  reassuring  that  there  is  a 
very  pronounced  two-zone  environment  based  on  temperature  measurement. 

A  hot  zone  up  at  the  ceiling  is  two  to  three  feet,  thick.  The  hot  gases 
are  recorded  by  the  thermocouples  at  8,  7  and  6  feet.  Temperatures 
recorded  by  the  remaining  thermocouples  in  the  lower  portion  of  the 
cabin  during  the  first  three  minutes  deviated  very  little  from  the 
ambient  temperature.  When  the  fire  developed  rapidly  due  to  the  col¬ 
lapse  of  the  panels  which  caused  burning  of  the  remaining  seats,  the 
two  temperature  zones  were  no  longer  apparent.  However,  there  was 
still  a  large  difference  in  temperatures  between  the  floor  and  the 
ceiling.  The  temperature  at  one-foot  and  two-foo*c  levels  were  less 
than  200  degrees,  whereas  the  temperature  at  the  ceiling  approached 
1000  degrees.  This  difference  in  temperature  was  reflected  by  the 
damage  of  the  interior  materials.  The  materials  in  the  upper  cabin 
were  virtually  destroyed,  whereas  those  near  the  floor,  especially 
the  carpet,  were  practically  undamaged.  The  carpet,  except  near  the 
fire  door  opening,  showed  very  little  damage. 
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-133  Temperature  Profile  with  Cabin  Materials 
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TIME  (SECONDS) 

Figure  6 


Figure  7  shows  the  hazard  development.  Instruments  were  installed 
inside  the  airplane  for  measuring  various  hazardous  components,  such 
as  smoke,  CC,  HCN,  temperature  and  0^  depletion.  Smoke  increased  early 
into  the  test;  the  oxygen  depletion  was  small  because  it  was  a  venti¬ 
lated  cabin.  The  predominant  hazards  were  CO  and  high  temperature. 

At  three  minutes,  the  CO  level  was  at  3000  parts  per  million,  HCN 
level  was  relatively  low  at  10  parts  per  million,  and  the  gas  tempera¬ 
ture  was  about  250°F  and  rising. 

This  concludes  ray  presentation.  I  hope  I  have  been  able  to  con¬ 
vey  to  you  the  major  characteristics  of  a  postcrash  aircraft  cabin 
fire  and  some  of  the  contrasting  features  with  residential  fires  which 
you  are  most  familiar  with. 

QUESTION: 

L&n  Cooper 3  National  Bureau  of  Standards.  You  have  very  much  down¬ 
played  the  pool  fire  aspect  of  a  hazard  and  I  wonder  if  you  could 
clarify  that  a  little  bit.  During  all  this  time,  you  saw  the  pool 
fire  going.  Earlier  you  showed  that  pool  fire  was  a  very  great 
hazard  in  and  of  itself.  Why  is  it  downplayed  in  this  scenario  and 
what  would  the  pool  fire  have  done? 

GUS  SARKOS: 

The  ultimate  goal  of  our  program  is  improved  test  methods  for  cabin 
materials.  Therefore,  we  are  just  trying  to  develop  a  realistic 
fire  scenario  which  uses  a  pool  fire  but  allows  the  interior  materials 
to  be  the  predominant  factor.  We  are  forcing  the  interior  materials 
to  be  the  predominate  factor  in  hazard  development  because  that  is  what 
we  are  interested  in.  We  do  not  want  to  mask  the  results  of  the  hazards 
developed  by  the  interior  materials  by  the  fuel  fire  hazards. 

QUESTION: 

Len  Coopert  National  Bureau  of  Standards.  Why  do  you  believe  this  to 
be  such  a  significant  scenario? 

GUS  SARKOS: 

We  are  focusing  the  scenario  to  come  out  that  way  because  we  are 
interested  in  the  materials.  The  accident  record  provides  very 
meager  statistics  to  derive  patterns  in  aircraft  accidents.  It  is 
very  difficult  to  come  up  with  a  typical  fire  scenario.  You  probably 
could  not  define  a  typical  fire.  We  derived  this  particular  scenario 
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FIGURE  7.  C-133  HAZARD  LEVELS  WITH  CABIN  MATERIALS  TEST  12/12/79 


which  corroborated  to  a  great  degree  the  actual  postcrash  accident 
that  I  have  just  talked  about.  There  was  an  intact  fuselage  with 
door  openings  adjacent  to  the  fire.  I  am  not  sure  I  have  answered 
your  questions.  Perhaps  we  can  get  together  later  on. 

QUESTION: 

Charles  Troha,  Consultant.  I  think  it  is  a  scenario  which  happens 
in  a  real  aircraft  fire.  The  question  is  what  effect  does  the  material 
have  on  the  total  involvement?  In  other  words,  you  developed  the 
scenario  for  a  pool  fire  and  were  you  able  to  subtract  any  of  that 
affect  to  show  the  real  affect  of  the  material? 

GUS  SAkKOS : 

If  I  didn't  mention  it,  I  meant  to  say  that  under  that  particular 
test  condition,  there  would  have  been  virtually  zero,  if  any,  hazard 
at  all  from  the  fuel  fire.  When  you  have  a  zero  wind  case,  with  a 
large  fuel  fire  next  to  the  opening,  you  get  very  little  accumulation 
of  hazards  from  the  fuel  fire,  It  is  hazardous  only  when  you  have 
flams  penetration.  That  particular  test  was  a  zero  wind  test.  The 
fuel  fire  hazards  were  minimal .  The  only  hazard  through  that  door 
opening  was  significant  radiant  heat.  A  flame  licking  in  randomly 
would  ignite  the  seat  as  it  was  being  cooked,  but  there  were  virtually 
no  hazards  from  the  fuel  fire  in  that  particular  test.  We  designed 
it  that  way.  We  did  not  want  to  mask  the  fuel  fire  hazards  from  those 
of  the  interior  materials. 
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FAA  MODELING  EFFORTS 


Thor  Eklund 

Fire  Safety  Branch 

FAA  Technical  Center 

First  of  all,  I  hope  the  weather  turns  out  mild.  For  many  of 
you  it  will  be  a  long  week  here.  Bob  Levine  and  Oliver  Foo  felt  that 
this  would  be  a  good  way  to  combine  a  number  of  different  efforts  and 
would  give  people  an  opportunity  to  see  this  area  and  cover  a  number 
of  topics.  There  will  be  ad  hoc  fire  mathematical  modeling  meetings. 
We  are  also  committed  to  have  a  workshop  on  the  DACFIR  model  developed 
by  Charles  MacArthur  under  FAA  contracts.  Because  our  efforts  in 
this  area  blossomed  over  the  last  year,  we  wanted  very  much  to  bring 
in  the  people  who  will  be  working  under  FAA  sponsorship.  It  will 
give  them  an  opportunity  to  learn  the  previous  work,  to  know  what  a 
postcrash  aircraft  fire  scenario  is,  and  to  distinguish  an  aircraft 
fire  from  a  home  or  dwelling  or  even  a  corridor  fire. 

In  1973  the  UDRI  contract  was  started.  This  was  the  FAA  pioneer 
effort  in  this  fire  modeling  area.  Over  the  years,  the  University  of 
Dayton  Research  Institute  (UDRI)  people  worked  closely  with  Boeing, 
NASA-Houston,  and  us  here  at  the  FAA  Technical  Center.  Until  the 
fall  of  1979,  this  was  under  sponsorship  of  the  Systems  Research  and 
Development  Services  (SRDS)  in  FAA  Headquarters  in  Washington.  The 
project  was  moved  under  our  sponsorship  in  1979.  We  felt  that  the 
SRDS  had  been  running  a  very  good  project  and  we  wanted  to  continue 
their  philosophy.  It  was  really  at  the  suggestion  of  Chuck  Troha 
that  we  started  an  interagency  agreement  with  the  National  Bureau  of 
Standards  (NBS)  to  further  go  into  the  field  modeling  as  well  as  zone 
modeling,  We  also  suggested  redirecting  the  Dayton  work  more  to  a 
postcrash  fire  scenario,  and  Dick  Kirsch  requested  our  involvement 
in  material  burning  at  Jet  Propulsion  Laboratory  (JPL)  under  Dr.  Kumar 
Ramohalli.  FAA  math  modeling  is  a  continuous  program  going  back  to 
1974.  Our  philosophy  on  modeling  is  very  much  the  same.  It  just 
happens  that  as  we  learn  more,  we  do  things  somewhat  differently. 


i'i 

There  are  five  tasks  that  concern  us  this  year.  First  of  all, 
we  are  interested  in  detailed  input  into  a  zone  model  like  DACFIR. 

A  zone  model  is  like  a  bathtub  upside  down  filling  up.  The  aircraft 
cabin  is  like  a  long  pipe  and  the  hot  layer  spreads  along  the  ceiling. 
There  are  considerable  heat  losses  and  changes  as  the  hot  layer  from 
a  fire  moves  along  the  cabin.  We  are  interested  in  further  develop¬ 
ment  of  smoke  layer  motion  at:  the  ceiling,  and  the  gas  dynamics  in 
detail . 


The  second  point  we  are  very  much  interested  in  is  thermal  impact 
at  openings.  This  has  been  fairly  well  described  in  what  you  saw 
earlier.  The  major  threat  from  the  iixternal  pool  fire  to  an  intact 
fuselage  is  through  an  opening,  The  equivalent  surface  temperature 
of  fire  is  of  the  order  of  1800  or  1900  degrees  F  and  the  radiative 
heat  flux  to  any  materials  inside  is  spectacular.  We  want  ultimately 
to  get  more  understanding  on  what  possibilities  we  would  have  of 
hardening  the  doorway.  We  also  want  to  know  at  what  rate  the  material 
fire  is  developed  in  that  area  from  such  huge  heat  fluxes. 


The  third  point  we  are  very  much  interested  in,  based  on  experi¬ 
mental  work  here  at  the  Technical  Center  over  the  last  three  years, 
is  the  effect  of  wind  on  fire  plume.  Given  an  external  pool  fire, 
the  wind  and  door  opening  configuration  is  the  predominant  factor 
regardless  of  the  material  involved.  We  have  asked  NBS  to  look  into 
the  pressure  distribution  around  the  fuselage  next  to  a  pool  fire. 

If  the  wind  is  blowing  over  a  fuselage  with  one  door  behind  the  fire 
and  the  other  door  facing  the  wind,  the  wind  will  drive  through  the 
aircraft  and  blow  everything  back  out  into  the  fire.  We  would  like 
to  get  some  quantitative  analyses  on  this  phenomenon.  Clearly,  when 
a  fire  burns  at  a  door,  the  stagnation  point  there  will  be  lost. 

We  don't  have  any  idea  of  the  magnitude  or  why  this  is,  but  this 
controls  the  ventilation  within  the  aircraft  during  the  fire.  We 
feel  that  this  is  a  very  important  point. 
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Item  number  four  that  we  are  working  closely  with  NBS  on  is 
correlation.  We  do  a  lot  of  small-scale  testing  and  we  do  a  lot  of 
large-scale  testing.  Like  everybody  else,  we  have  problems  with  cor¬ 
relating  small-scale  tests  to  large-scale  tests.  This  is  on  our  mind 
very  much  right  now. 

The  fifth  item  that  we  are  looking  at  and  are  very  interested  in 
is  actually  the  mechanisms  of  burning  and  flame  spreading  over  air¬ 
craft  type  materials.  We  are  working  through  NASA  and  with  NBS  on 
this  topic.  It  is  hard  enough  to  study  the  flame  spread  and  burning 
mechanism  with  a  simple  and  uniform  material  like  plexiglass.  The 
aircraft  materials  have  fire  retardants  and  are  often  laminated  one 
way  or  another.  The  existing  data  and  test  methods  are  insufficient 
for  aircraft  materials.  It  is  mandatory  that  we  get  some  answers 
soon  in  this  area. 

Those  are  the  five  areas  that  we  are  interested  in  at  this  time 
in  math  modeling,  I  would  like  to  reinforce  some  things  that  Gus 
Sarkos  presented  to  you.  We  are  spending  something  on  the  order  of 
15  to  20  percent  of  cur  cabin  fire  safety  budget  on  math  modeling 
because  it  is  important.  We  have  a  lot  of  other  high  priority 
obligations.  An  aircraft  cabin  has  the  shape  of  a  long  tube  and  is 
packed  with  plastic  materials  (side  walls,  ceiling,  carpet,  and  seats). 
Furthermore,  it  is  densely  populated  with  people.  A  huge  heat  source 
is  an  external  pool  fire.  I  will  give  you  an  example  of  how  severe 
this  pool  fire  can  be.  When  we  did  small-scale  tests,  we  ignited  a 
pan  with  around  five  gallons  of  fuel,  I  could  not  get  close  to  that 
pan.  A  747  taking  off  with  its  total  fuel  has  close  to  50,000  gallons. 
I  can't  in  ir.y  mind  imagine  what  kind  of  fire  you  could  make  with  that 
amount  of  fuel.  That  is  something  we  want  to  hammer  in.  It  is  a  very 
serious  heat  source.  We  do  know  wind  and  door  openings  are  important 
but  we  are  not  too  comfortable  with  our  understanding  of  the  relation¬ 
ships  to  fire. 
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There  are  two  different  aircraft  fires.  The  in-flight  fire  problem 
is  to  put  it  out  before  it  gets  severe.  But  if  it  gets  severe,  nobody 
can  get  out.  The  postcrash  fire  is  a  very  rapid  developing  fire  and 
everybody  has  to  get  out  very  quickly. 

I  would  also  like  to  say  a  few  things  about  the  aircraft  materials 
burning  phenomena.  The  composite  panels  do  not  burn  so  well.  They  may 
cook  out  and  disintegrate  and  fall.  Other  than  the  panel  outer  layer, 
the  rest  of  the  panel  components  do  not  burn  heavily.  The  carpet,  as 
long  as  the  fire  is  not  coming  from  underneath  the  aircraft  or  from 
the  cargo  area,  is  generally  pretty  flame  resistant.  The  urethane 
seats  with  various  coverings  right  now  seem  to  be  the  big  factor  in  a 
fire.  The  materials  which  are  placed  in  different  orientations  suffer 
from  different  exposure  conditions.  Their  behaviors  are  quite  differ¬ 
ent  even  under  the  same  overall  test  conditions. 

There  are  two  final  remarks  about  the  aircraft  configuration.  The 
cabin  is  a  longitudinal  ore.  Great  buoyancy  forces  are  difficult  to 
be  generated,  as  compared  to  that  from  an  enclosure  fire.  Also,  it  is 
pressurized  for  the  in-flight  fire.  These  can  have  effects  on  an  analysis. 

I  would  like  to  say  a  few  words  now  about  what  we  are  ultimately 
looking  for.  Like  everybody  else,  we  would  like  to  have  perfect  math 
models  which  would  predict  everything.  We  really  don't  believe  that 
will  be  the  case.  Right  now  we  are  interested  in  separately  looking  at 
gas  dynamic  development  and  material  burning  phenomena.  In  the  future, 
we  would  like  those  to  be  bridged.  That  is,  once  you  know  what  is 
happening  within  the  fuselage,  you  can  start  saying  what  you  know  about 
materials  behaviors  under  various  exposures.  That  is  a  little  further 
downstream. 

We  are  interested  now  in  whatever  test  models  are  available  and 
use  them  in  our  test  programs.  We  want  to  know  that  we  are  making  the 
right  measurements  in  the  right  places.  The  FAA  wants  to  develop  new 
standards,  which  have  to  be  very  defensible.  Our  C-133  test  represents 
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only  one  scenario.  One  way  or  another  the  validity  of  that  scenario 
has  to  be  demonstrated  either  by  mathematical  solutions/small-scale 
tests  or  by  logical  arguments.  We  need  mathematical  modeling  to  give 
us  a  better  handle  on  the  type  of  situations  that  can  occur  and  possibly 
what  we  can  do  about  them. 

We  are  interested  in  mechanical  type  countermeasures.  For  instance, 
in  the  case  of  the  wind  caused  pressure  distribution,  we  could  have 
ventilation  countermeasures.  In  the  case  of  better  elucidation  of  fire 
in  the  doorway,  we  might  be  able  to  give  a  rational  basis  for  fire 
hardening  procedures.  If  we  can  get  a  better  handle  on  the  mechanisms 
of  burning,  we  might  be  able  to  design  better  material.  There  is  one 
effort  we  are  involved  in  now  that  might  lead  to  such  a  solution. 

In  summary,  we  need  the  modeling  to  expand  our  scenarios,  to  find 
the  key  test  parameters  that  we  should  be  looking  at,  and  to  correlate 
small-  and  full-scale  tests.  I  would  like  to  say  aircraft  materials 
are  very  good  now.  It  is  the  magnitude  of  the  postcrash  fuel  fire  and 
also  the  lack  of  egress  capability  in  an  in-flight  situation  that  makes 
the  aircraft  fire  still  a  terrifying  situation  to  think  about.  Any 
kind  of  elucidation  we  can  get  theoretically  or  experimentally  I  hope 
we  can  put  to  good  use. 
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The  DACFIR  Model,  Dayton  Aircraft  Fire  Model,  is  in  its  third 
version.  Even  though  this  program  has  been  going  on  for  some  time, 
this  third  and  final  version  was  just  developed  within  the  last  three 
months.  It  still  has  some  testing  and  perhaps  debugging  to  be  done, 
but  the  results  to  date  are  very  encouraging.  The  handout  packages 
contain  equations,  assumptions,  results,  etc,,  all  pertaining  to  this 
third  version. 

The  computer  program  at  this  point  is  in  good  condition  and  can 
be  distributed  for  those  who  wish  to  get  a  copy  of  it.  I  believe 
that  FAA  will  make  these  available  in  a  very  short  time.  The  report 
on  this  third  version  of  the  program  and  the  computer  listing,  in  a 
tabular  form,  will  be  available  within  a  month  or  a  month  and  one- 
half  after  the  FAA  review  is  completed. 

At  the  very  start  of  this  program,  the  specifications  for  a  com¬ 
puter  model  on  aircraft  cabin  fires  were  laid  out  in  the  statement  of 
work.  The  objective  of  the  model  (Figure  1)  was  to  assess  the  smoke 
and  toxic  gas  accumulation  in  the  cabin  resulting  from  an  exterior 
fire.  As  you  have  heard  earlier,  the  situation  has  changed.  The  FAA 
is  more  interested  in  the  exterior  fire  and  its  effect  on  the  inter¬ 
ior  materials.  When  we  started  the  program,  the  emissions  scenario 
was  an  interior  ignition  which  might  be  a  ruptured  fuel  line  through 
the  floor  or  a  spilled  flammable  liquid  in  the  interior,  and  the 
effect  of  the  interior  material  on  survivability.  We  did  not  formulate 
the  problems  by  starting  with  the  first  principles  of  thermodynamics. 

It  was  not  possible  then  and  still  may  not  be  possible  now,  We  were 
looking  for  a  practical  first-cut  engineering  solution  to  predicting 
the  survivability  of  the  cabin.  The  emphasis  was  on  the  practical 
method  and  the  method  that  could  be  used  for  safety  decision  making. 


OBJECTIVE  OF  THE  CABIN  FIRE  MODEL 


A  very  important  part  of  our  model  is  that  it  uses  the  materials 
flammability  behavior  from  the  laboratory  scale  test  data.  It  was 
not  clear  in  the  objective  of  the  model  which  particular  test  result 
should  be  used  so  we  had  some  leeway  in  choosing  the  test  result. 

The  model  was  a  way  of  correlating  mathematically  the  small-scale 
test  with  the  large-scale  test. 

Figure  2  shows  a  sketchy  ove  view  of  what  the  computer  model 
does.  The  model  takes  input  information  about  the  cabin,  about  the 
materials,  and  about  the  particular  ignition  scenario  as  input  and 
produces  histories  of  the  gas  composit ions  and  temperatures  of  the 
cabin  atmosphere.  It  does  predict  in  a  simple  fashion  the  regions 
of  fire  spread  damage.  The  model  consists  of  two  main  parts.  A  fire 
spread  part  which  attempts  to  predict  the  amounts  of  materials  that 
are  affected  by  the  fire  burning  as  a  function  of  time,  i.e.,  the 
development  of  the  interior  fire.  The  other  main  part  is  gas  dynam¬ 
ics  simulating  the  cabin  atmosphere  motion  and  combustion  products 
accumulation  in  the  atmosphere  as  the  fire  progresses. 

cabin  atmosphere  model  in  version  three  is  a  zone  model.  As 
shown  in  Figure  3,  the  cabin  is  divided  into  two  zones.  The  upper 
zone  contains  the  accumulation  of  the  combustion  products.  The  lower 
zone  is  primarily  uncontaminated  air  which  gets  pumped  into  the 
upper  through  the  actions  of  flames  and  plumes,  etc.,  in  the 

cabin.  The  lower  picture  in  Figure  3  shows  a  rough  cross  section  of 
the  cabin  to  emphasize  that  we  compute  temperatures  and  species  mass 
fractions.  We  do  computations  on  the  temperature  and  composition  for 
both  zones.  This  is  not  really  that  important  for  the  lower  zone  in 
this  version,  but  the  program  is  set  up  to  handle  any  future  develop¬ 
ment  that  might  occur,  in  particular,  mixing  between  the  zones. 

The  fire  spread  part  of  the  model,  covers  a  section  of  the  cabin. 
We  do  not  attempt  to  predict  the  involvement  of  the  materials  over 
all  of  the  interior  surfaces.  This  is  mainly  to  economize  the  computer 
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DAYTON  AIRCRAFT  CABIN  FIRE  MODEL 
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code  and  run  time.  Also,  in  most  of  the  situations  when  the  interior 
materials  are  involved,  the  flame  will  not  spread  over  the  materials 
beyond  two  or  three  seat  row  sections  from  the  origin  of  the  fire 
before  the  interior  of  the  cabin  is  really  not  inhabitable  at  all. 

We  are  not  trying  to  predict  the  development  of  the  fire  up  to  flash- 
over  or  a  fully  developed  compartment  fire,  but  only  that  first  three, 
five  or  maybe  ten  minutes  during  which  people  may  still  be  able  to 
escape  and  conditions  haven't  become  intolerable  yet.  It  is  the 
objective  of  the  model  to  predict  the  time  at  which  the  cabin  be¬ 
comes  unsurvivable. 


Figure  4  shows  a  very  crude,  but  effective,  presentation  of  the 
cabin  interior.  The  seats  consist  of  horizontal  and  vertical  planes 
in  L-shape,  The  surface  of  these  planes  is  divided  into  square  regions 
which  are  named  fuel  elements.  The  method  of  computing  fire  spread 
in  the  program  is  through  a  method  of  tracking  these  elements  from  the 
undisturbed  state  into  a  flaming  state  or  a  smoldering  state  then  into 
a  burned  out  state,  etc.  It  is  a  discrete  step-byrstep  description 
of  fire  behavior  from  the  materials.  It  is  an  oversimplification, 
but  a  good  first  cut,  in  handling  the  very  complex  geometry  of  a 
cabin  interior  with  furnishings. 


The  cabin  interior  surfaces  were  divided  into  square  regions. 

The  dimension  of  these  squares  was  one-half  foot,  mainly  because  it  was 
a  convenient  length  scale  for  the  interior  of  the  cabin  and  also  it 
did  not  really  create  an  excessive  amount  of  computer  storage.  The 
program  could  be  refined  to  have  smaller  element  sizes  to  predict  areas 
more  precisely,  but  one-half  foot  is  a  good  practical  compromise  right 
now,  < 

Figure  5  shows  how  the  development  of  the  fire  is  tracked  in  the 
computer  model  by  adding  the  shaded  squares  which  are  regarded  as  being 
on  fire  and  being  a  source  of  heat  and  smoke  and  gas  emissions.  These 
elements  (burning  or  burned  out)  are  determined  by  the  particular 
material  data  supplied  to  the  program  inputs.  These  geometric  regions 


52 


TRACKING  THE  FIRE  SPREAD 


Groups  of  Flaming  Elements  form  Fire  Bases 


then  are  the  so-irces  of  combustion  that  is  fed  into  a  subprogram  cabin 
atmosphere  to  determine  smoke  development  and  gas  compositions. 

DACFIR  '’orsion  3  has  a  major  refinement  and  improvement  over  the 
earlier  versions.  As  shown  in  Figure  6,  the  model  is  designed  to 
simulate  fire  in  a  cabin  which  has  one  to  four  compartments.  Earlier 
versions  just  considered  the  cabin  as  one  long  room.  In  DACFIR-3, 
one  room  can  be  divided  into  four  rooms  attached  linearly  along  the 
cabin  and  each  of  these  compartments  may  have  one  to  six  vents  or  doors 
or  escape  hatches  or  openings  to  the  exterior  or  through  the  dividers 
to  one  another.  DACFIR-3  has  retained  the  capability  in  DACFIR-2  of 
handling  prescribed  flows  at  doors.  This  was  used  to  compare  the  model 
to  test  data  in  which  one  or  more  of  the  doors  had  a  forced  flow  from 
the  floor,  The  user  is  allowed  to  specify  the  temperature  and  composi¬ 
tion  of  this  inflow  gas  (at  least  one  vent) ,  This  is  a  first  step  in 
being  able  to  have  the  model  simulate  the  effect  of  the  exterior  fire. 

The  computer  program  in  DACFIR-3  is  considerably  different  from 
the  earlier  versions  (Figure  7) .  The  earlier  versions  used  a  very 
primitive  method  of  integrating  the  equations  of  the  model  cabin 
atmosphere.  A  very  good  technique  which  was  used  in  the  Harvard 
Computer  Fire  Code-3  is  implicit  (trapezoidal  rule)  integration  of 
the  atmosphere  equation  using  the  Newton-Raphson  technique,  I  am 
really  surprised  at  the  stability  and  reasonable  economy  of  this  tech¬ 
nique  over  the  other  integration  methods.  DACFIR-3  adopted  a  modular 
construction,  at  least  in  the  cabin  atmosphere  part  of  the  program. 

The  model  was  designed  along  the  lines  of  the  Harvard  code  and  other 
codes  developed  in  the  fire  mathematical  modeling  workshop  group. 

The  subroutines  that  contain  the  modules  can  be  independently  re¬ 
moved  and  replaced  if  necessary.  DACFIR-3  is  a  computer  program 
that  is  easy  to  maintain,  and  will  be  easy  to  upgrade  when  future 
improvements  in  zone  modeling  are  available. 

An  overall  flow  chart  of  the  computer  code  is  shown  in  Figure  8. 
Essentially,  there  are  two  parts,  i.e.,  the  flame  spread  part  and  the 
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DACFIR  VERSION  3 
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Figure  8 


gas  dynamics  part.  Subroutine  "ATMOS"  is  the  controlling  part  of 
the  computer  program  for  the  cabin  atmosphere  model.  The  differential 
equations  of  fire  physics  were  integrated  at  a  relatively  small  time 
step.  A  cumulative  structure  for  the  program  was  used.  The  integra¬ 
tion  of  fire  physics  equations  can  be  advanced  independently  relative 
to  the  flame  spread  subroutines.  The  flame  spread  subroutines  scan 
all  elements  and  determine  which  elements  are  ignited  and  which  ele¬ 
ments  are  burned  out,  etc. 

DACFIR-3  uses  a  zone  model  approach  which  is  patterned  after  the 
model  developed  by  Prof.  Emmons  and  Dr.  Mitler  of  Harvard  and  Dr. 
Quintiere's  model  at  NBS.  DACFIR-3  uses  a  two-zone  concept  to  model 
fire  enclosure.  It  deals  with  multiple  compartments.  Figure  9  shows 
two  compartments  in  a  cabin.  Each  compartment  has  an  upper  and  a 
lower  zone.  The  computer  program  will  allow  fires  to  exist  within  the 
lower  zone  or  the  upper  zone.  It  is  convenient  that  the  upper  zone 
was  brought  down  through  the  plume  to  the  fire  base.  Dr.  Quintiere 
documented  this  idea.  From  a  conceptual  standpoint,  it  minimizes  the 
problem  about  the  interface  between  the  plume  and  the  zone.  The 
variables  for  these  gas  zones  are  temperatures,  density,  and  the  composi¬ 
tions.  Shown  in  Figure  9  is  a  particular  case  where  a  flow  exists, 
not  only  between  compartments  on  the  right  but  also  on  the  left  through 
an  open  door  to  the  exterior,  A  flow  out  of  the  door  and  return  flow 
into  the  lower  zone  are  also  seen  in  the  picture. 

Figure  10  is  a  list  of  the  variables  of  the  atmospheric  model. 

In  particular,  a  pressure  for  the  entire  compartment  is  calculated 
and  used  as  a  single  reference  pressure.  Conservation  equations  are 
shown  in  Figure  11. 

Another  new  feature  of  the  DACFIR  model  is  the  very  simple  global 
one-step  model  of  combustion  chemistry  which  is  shown  in  Figure  12. 

One  of  the  problems  in  testing  is  an  understanding  of  what  the  mass 
fraction  of  water  vapor  might  be  in  the  gas.  In  certain  situations. 
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Gas  Zone  Control  Volumes  and  Mass  Flows 


Figure 


variari.es  of  the  cabin  atmosphere  model 


Variable 

Symbol 

Lower  zone  species 
mass  fractions 
(j  values)* 

xij 

Upper  zone  species 
mass  fractions 
(J  values) 

xuJ 

Pressure 

Pf 

Lower  zone  density 

°* 

Upper  zone  density 

pi 

Lower  zone  temperature 

t[ 

Upper  zone  temperature 

Ti 

Lower  zone  volume 

V1 

Upper  zone  volume 

V1 

u 

Thermal  discontinuity 

71 

position 

<-d 

Materials  surface 
temperature** 

T  1 

rsk 

*  Minimum  value  of  j  is  5  and  the  maximum  11 
##  Minimum  value  of  k  is  1  and  the  maximum  20  per  compartment 


Figure  10 
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CONSERVATION  EQUATIONS 
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Figure  11 
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OXYGEN  CONSUMPTION  AND  PRODUCT  EMISSION 
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Figure  12 


some  acid  gases  might  be  scrubbed  out  by  condensing  water  vapor. 

This  inspired  the  idea  of  using  a  very  simple  combustion  model  for 
all  fuels,  even  the  polymers  that  contain  constituents  other  than 
carbon  and  hydrogen.  This  reaction  is  the  source  term  in  the  species 
conservation  equations  for  oxygen,  nitrogen,  CO^,  and  H^O,  and  fuel 
vapor . 

DACFIR-3  deals  with  only  four  species  with  the  fuel  vapor  mass 
fraction  set  to  zero.  It  is  not  assumed  that  there  is  any  unburned 
fuel  vapor  existing  in  the  plume  or  in  the  upper  and  lower  zones. 

It  is  assumed  that  immediately  as  the  fuel  vapor  touches  the  surface 
of  the  upper  zone  control  line,  it  is  completely  reacted  and  the 
products  are  carried  throughout  the  upper  zone.  The  computer  program 
is  structured  to  have  a  non-zero  fuel  vapor  mass  fraction  in  the  upper 
zone. 

One  of  the  unfortunate  things  we  don't  know  about  is  a  measure  of 
mass  burning  rate  as  a  function  of  anything.  The  test  data  that  we 
use  is  a  derivative  of  that  fundamental  quantity  in  terms  of  heat 
release  rate  and  product  release  rate.  We  can  estimate  what  that 
mass  burning  rate  is  and  use  that  in  the  species  terms  by  taking  the 
heat  release  rate  and  dividing  it  by  heat  combustion  for  materials. 

Flame  and  plume  entrainment  is  a  problem  that  has  haunted  zone 
modelers  for  a  long  time.  At  a  lower  zone,  air  is  entrained  into  the 
plume.  It  travels  in  the  upper  zone  and  dilutes  the  combustion  pro-, 
ducts.  An  air  entrainment  model  first  introduced  by  Prof.  Steward  in 
Combustion  Science  and  Technology.  1970,  and  refined  by  Dr.  Fang  of 
NBS  is  used  in  DACFIR-3,  shown  in  Figure  13.  This  model  does  differ¬ 
entiate  between  the  combusting  zone  with  heat  generation  and  the 
plume  without  heat  generation  above  combustion  zone.  There  are  two 
entrainment  constraints.  The  mathematical  formulations  are  also 
given  in  Figure  13.  This  is  a  classic  example  of  a  part  of  the  model 
which  can  be  removed  easily  as  one  subroutine  and  could  be  replaced 
with  another. 
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FLAME  AND  PLUME  ENTRAINMENT 


Figure  13 


The  other  major  new  item  in  DACFIR-3  is  a  vastly  improved  method 
of  calculating  the  pressure  and  buoyancy  driven  flows  through  com¬ 
partment  vents.  This  is  a  method  that  has  been  used  by  others,  par¬ 
ticularly  Harvard  researchers  and  Prof.  Tanaka.  The  basic  idea  is  to 
use  the  hydrostatic  law  to  compute  the  pressure  variation  across  from 
either  side  of  the  vent  opening  and  then  compute  the  flow  through  the 
opening  from  the  orifice  equation  knowing  the  pressure  difference. 

In  earlier  zone  models,  one  fixed  the  pressure  in  the  floor  to  be  the 
ambient  pressure  and  all  of  the  flow  or  most  of  the  plume  was  due  to 
the  difference  of  pressures  due  to  buoyancy.  In  DACFIR-3,  the  pres¬ 
sure  at  the  floor  is  a  variable  also.  A  pressure  difference  between 
two  compartments  at  the  floor  could  generate  a  pressure-driven  flow. 
When  there  is  a  density  difference  between  two  zones  and  the-  rate  of 
change  of  pressure  with  height  differs,  it  leads  to  a  situation  where 
a  neutral  plane  exists  in  the  doors. 

Figure  14,  (a)  and  (b) ,  shows  flow  from  a  single  compartment  to 
the  exterior  while  (c)  and  (d)  show  jws  be  ween  two  compartments 
with  the  pressure  profiles  intersecting  at  veral  points.  It  means 
that  there  is  a  flow  from  compartment  one  to  compartment  two  above  a 
height  and  a  reverse  flow  below.  It  is  possible  to  have  two  flows 
between  the  upper  zones,  This  is  possible,  but  I  can't  say  I  have 
ever  seen  it  in  any  of  our  test  runs.  I  am  not  sure  anyone  ever  will, 
but  the  program  is  set  up  to  handle  this  very  complicated  situation. 

Figure  15  shows  vent  flow  computations.  The  formulation  is  no 
different  from  that  presented  by  Prof.  Emmons  in  1978.  Our  method  of 
solution  is  a  little  unique.  Rather  than  a  very  complicated  logic 
tree  to  select  certain  formulas,  we  have  taken  the  hydrostatic  law 
which  is  the  pressure  as  a  linear  function  of  height,  that  breaks  at 
the  thermal  discontinuity  position  in  each  compartment.  Take  P  as  a 
function  of  Z  and  find  the  intersections  of  the  pressure  profiles  in 
each  compartment.  We  solve  those  equations  and  then  we  decide  whether 
those  neutral  planes  are  physically  possible.  For  example,  some 


PRESSURE  AND  BUOYANCY  DRIVEN  FLOWS  THROUGH  COMPARTMENT  VENTS 
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Figure  14 


intersections  may  occur  below  the  floor  level  or  may  occur  at  very 
high  values  of  Z--higher  than  the  ceiling  of  the  compartment.  Once 
we  have  decided  whether  there  are  any  neutral  planes  in  the  door  then 
we  also  check  for  the  position  of  the  thermal  discontinuity.  The  re¬ 
sultant  pressure  difference  is  a  piecewise  linear  function  of  the 
heights  across  the  vent,  either  between  two  compartments  or  from  one 
compartment  to  the  exterior.  P  is  a  function  of  Z  and  this  function 
is  either  constant  when  the  densities  are  the  same  on  either  side 
or  is  a  linear  function  of  Z.  In  either  case,  it  can  be  integrated 
analytically  to  obtain  two  expressions  for  the  mass  flow  rate 
which  is  as  a  function  of  the  pressure  difference.  In  the  tests 
we  have  run  so  far,  it  ran  pretty  well.  The  whole  subroutine  was 
programmed  into  100  lines  of  PORTRAN  code. 

QUESTION: 

Dr.  Michael  A.  Delichatsios,  Factory  Mutual  Research  Corporation. 

How  do  you  know  the  pressure  distribution? 

CHARLES  MacARTHUR: 

I  know  the  pressure  at  the  floor,  and  I  know  the  density  and  so  all 
I  need  is  the  hydrostatic  law  to  predict  the  pressure. 

QUESTION: 

Dr.  Michael  A.  Delichatsios ,  Factory  Mutual  Research  Corporation. 

If  such  flows  exist,  there  are  a  lot  of  eddies.  Would  this  change 
the  flow  completely? 

CHARLES  MacARTHUR: 

Yes. 

QUESTION: 

Dr.  Michael  A.  Delichatsios ,  Factory  Mutual  Research  Corporation. 
Would  the  hydrostatic  flow  equations  apply? 

CHARLES  MacARTHUR: 

No.  The  hydrostatic  law  does  not  apply  in  the  case  where  there  is 
any  velocity  at  all  in  reality.  This  is  an  approximation  that  needs 
relatively  low  speed  flows.  We  can  approximate  the  true  pressure 
distribution  by  the  hydrostatic  law. 
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QUESTION: 

Dr.  Michael  A,  Deiiahatsios  *  Faatory  Mutual  Research  Corporation. 

Can  a  zone  model  take  into  account  something  that  is  really  the 
providence  of  a  field  model--that  is  a  precise  calculation  of  the 
nressure  and  velocity  fields  at  each  point. 

CHARLES  Mac  ARTHUR: 

The  zone  models  can’t  do  that,  but  you  brought  up  an  important  fact 
that  I  might  have  missed.  There  is  work  now  being  done  at  NBS  on 
the  iiixing  at  the  door  where  two  relatively  high  speed  flows  in 
opposite  directions  with  a  very  high  shear  rate  occur  at  the  thermal 
end  of  the  neutral  plane.  There  are  eddies  which  promote  mixing. 

That  is  how  some  of  the  upper  zone  product s--temperature  and  species 
get  mixed  aown  in  the  lower  zone.  DACFIR-3  does  not  have  a  model 
of  that  mixing  because,  as  T  understand,  there  isn't  really  a  model 
available.  One  of  the  reasons  we  do  mass  species  and  energy  balances 
on  the  lower  zone  is  anticipation  of  having  a  mixing  at  the  door 
parameterized  in  a  formula  available  so  we  can  predict  species  con¬ 
centrations  in  the  lower  zone. 

QUESTION: 

Dr.  John  de  Ris ,  Faatory  Mutual  Research  Corporation.  I  would  just 
like  to  comment.  Mike  brings  up  a  point  that  I  suspect  is  not  an 
issue,  but  I  think  it  has  to  get  settled  by  looking  at  the  Richardson 
number  of  those  flows.  One  may  be  able  to  estimate  at  least  in  some 
crude  way  what  the  Richardson  numbers  would  be  here  and  I  wonder  if 
anyone  has  done  that? 

CHARLES  MacARTHUR: 

We  knew  this  was  in  the  wind  and  that  is  why  we  structured  the  third 
version  this  way.  This  will  be  discussed  at  the  workshop  on  Friday 
and  we  will  give  you  experimental  data. 

Our  models  of  the  heat  transfer  from  the  upper  zone  are  admit¬ 
tedly  very  simplified,  shown  in  Figure  16,  because  we  can't  spend  a 
lot  of  computing  time  and  effort  into  developing  individual  parts. 

The  convective  flow  of  the  heat  zone  is  just  a  simple  constant  film 
coefficient  multiplying  the  area  of  the  surface  to  which  the  convec¬ 
tion  is  taking  place  and  difference  in  temperature.  I  think  the 
assumed  sign  is  incorrect  in  Figure  16,  The  gas  temperature  would  be 
higher  than  surface  temperature  in  most  situations,  at  least  in  the 
upper  zone,  so  Q  would  be  negative.  In  the  computer  code,  the  signs 
are  kept  in  the  right  fashion.  A  major  refinement,  even  though 
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CONVECTION  AND  RADIATION 
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Figure  16 


convective  loss  to  the  surfaces  is  not  one  of  the  major  heat  loss 
terms  in  the  energy  equation,  would  be  to  have  some  better  estimate 
of  what  the  convection  coefficient  is.  We  have  adopted  a  typical 
value  for  turbulent  flow  over  a  flat  surface  used  particularly  by 
Harvard.  It  is  up  in  the  air  as  to  whether  it  is  any  good.  Fortu¬ 
nately,  the  model  is  not  too  sensitive  to  convective  heat  losses. 

If  we  want  to  do  some  very  careful  analysis  of  the  temperatures  of 
materials  and  use  the  temperature  of  surfaces  to  predict  flame 
spread  rate,  then  we  have  got  to  go  back  and  look  at  this  convective 
loss  terra  a  little  closer. 

Radiation  loss  is  one  of  our  larger  terms  in  most  scenarios. 

The  radiation  absorbed  by  the  gas  zone,  Qrin»  takes  into  account 
surfaces  lining  the  zone  in  the  first  term  and  radiation  from 
neighboring  zones  in  the  second  term.  Radiation  emitted  by  the  zone 
was  calculated  by  using  grey  gas  approximation  with  a  mean  beam  length 
approximation.  This  is  an  equation  which  first  appeared  in  Dr. 

Quint iere's  work  in  estimating  emittance  using  the  smoke  density  and 
also  an  absorption  coefficient  for  gas  species  and  gas  band  radiation. 
The  Qrout  term  which  is  the  total  radiated  energy  by  the  upper  zone 
to  everything  with  Asurf  being  the  total  upper  zone  surface  area. 

In  Figure  17,  the  equation  of  state  or  the  gas  law  for  each 
specie  is  given  in  terms  of  partial  pressure  or  density.  The  con¬ 
servation  of  volume  for  each  cell  and  the  interface  height  at  the  dis¬ 
continuity  are  also  given.  The  foregoing  physics  equations  are  used 
to  calculate  the  gas  dynamics  in  the  cabin. 

The  numerical  procedures  for  solving  these  equations  are  out¬ 
lined  in  Figure  18.  Trapezoidal  rule  integration  of  ordinary  dif¬ 
ferential  equations  are  coupled  with  the  Newton-Raphson  iterative 
method  for  a  set  of  algebraic  equations.  This  technique  follows 
the  latest  developments  by  the  Harvard  University  Fire  Research 
Group  and  is  very  successful  in  terms  of  numerical  stability  and 
computer  time  usage. 
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Figure  17 


NUMERIC  SOLUTION 

OF  THE  CABIN  ATMOSPHERE  EQUATIONS 
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Figure  18 


In  order  to  calculate  the  surface  temperatures  of  materials, 
three  assumptions  were  made,  as  shown  in  Figure  19.  The  material 
properties  are  assumed  homogeneous  and  the  surface  temperatures  are 
assumed  constant  during  an  integration  step.  A  simple  Euler  integra¬ 
tion  scheme  is  used  to  integrate  the  energy  balance  equation  at  the 
material  surface  (Figure  20). 

The  gas  temperature  becomes  high  and  gives  out  heat  to  materials 
ahead  and  in  the  lower  zone.  These  equations  are  given  in  Figure  21. 
The  view  factors  are  given  in  cabin  geometry  dimensions  which  are 
indicated  in  Figure  22.  A  Cartesian  coordinate  system  was  used  in 
the  model  for  convenience.  Three  dimensional  indices  are  used  to 
label  the  cells. 

The  seven  element  states  are  defined  in  Figure  23.  The  allow¬ 
able  transitions  from  one  element  state  to  another  element  state 
are  shown  in  Figure  24.  The  computer  code  has  a  subroutine  to  de¬ 
termine  which  element  state  each  cell  is  in.  An  element's  transi¬ 
tion  from  one  state  to  another  is  governed  by  the  properties  of  the 
material  associated  with  che  element  and  by  the  element's  relation¬ 
ship  to  the  fire  in  the  cabin.  A  fire  is  defined  by  a  set  of  con¬ 
tinuous  elements  in  state  3. 

The  rate  at  which  a  flame  front  propagates  depends  upon  several 
factors.  The  factors  considered  in  this  program  are  the  type  of 
material  at  the  edge  of  a  fire,  the  size  of  fire,  orientation  of  the 
surface,  and  the  background  radiation  level.  The  flame  spread  rates 
for  a  given  materia)  are  input  data  to  the  computer  program  and  are 
in  a  tabular  form  as  functions  of  heat  flux  (Figure  25) .  The  heat 
flux  to  elements  adjacent  to  flaming  elements  is  calculated  based  on 
the  size  of  an  adjacent  fire  and  the  overall  background  radiation 
level.  Three  flame  spread  rates  are  associated  with  a  vertical  sur¬ 
face:  vertical  up,  vertical  down,  and  horizontal.  One  flame  spread 
rate  is  associated  with  horizontal  surfaces.  The  rates  and  directions 
are  shown  in  Figure  26. 
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MATERIALS  SURFACE  TEMPERATURES 
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Figure  19 


MATERIALS  SURFACE  TEMPERATURES 
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Figure  20 


UPPER  ZONE  RADIATION  TO  MATERIALS 


Figure  21 


ELEMENT  STATES 


State  1  -  VIRGIN 

The  element  Is  in  its  virgin  state;  it  has  not 
been  directly  affected  by  the  fire. 

State  2  -  SMOLDERING 

The  element  is  undergoing  nonflaming 
decomposition. 

State  3  -  FLAMING 

The  element  is  undergoing  self-sustaining 
combustion. 

State  4  -  CHARRED 

The  element  has  burned  out  and  will  no  longer 
smolder  or  burn. 

State  5  -  HEATING,  NO  FLAME  CONTACT 

The  element  is  receiving  heat  flux  sufficient 
to  cause  it  to  smolder  but  smoldering  has  not 
yet  begun. 

State  6  -  HEATING,  WITH  FLAME  CONTACT 

The  element  is  being  touched  by  the  flames  of 
a  fire  but  has  not  yet  ignited. 

State  7  -  SMOLDERING  AND  COOLING 

The  element  began  smoldering  when  the  heat  flux 
reached  a  specified  level;  the  flux  has  now 
dropped  below  that  level  but  the  material  is 
still  smoldering. 


Figure  23 
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INPUT  DATA  ON  MATERIALS  FIRE  BEHAVIOR 
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Figure  25 


The  flame  radiation  calculation  is  based  on  an  equation  derived 
by  Dayan  and  Tien,  Combustion  Science  and  Technology,  Vol .  9,  1974 
pp.  41-47,  for  a  cylindrically  shaped  fire  on  a  horizontal  surface 
facing  upward.  It  is  used  in  the  present  model  to  compute  the 
radiation  level  at  the  edge  of  any  fire  base.  These  equations  and 
the  equations  used  for  calculating  flux  levels  are  shown  in  Figure  27. 
The  flame  height  is  calculated  with  the  equation  derived  by  Steward 
and  Fang  and  refined  by  Fang  (NBSIR  73-115).  The  smoldering  range 
is  obtained  from  the  model  by  Dayan  and  Tien,  shown  in  Figure  28. 

The  statistics  of  computer  program  for  DACFIR-3  are  given  in 
Figure  29.  There  are  4050  source  statements  with  a  required  memory 
core  of  325,000  bytes.  The  execution  time  on  a  DEC  VAX-11/780  is 
1500  seconds  CPU  time  for  a  simulated  time  of  400  seconds.  The  sample 
outputs  are  given  in  Figures  30-32.  The  test  cases  to  be  simulated 
are  three  test  runs  performed  in  a  737  fuselage  at  Johnson  Space 
Center/NASA.  The  test  conditions  are  described  in  Figure  33. 

The  height  of  the  thermal  discontinuity  is  given  in  Figure  34. 

As  time  goes  by,  the  thermal  discontinuity  descends  down  to  a  lower 
level  as  the  upper  layer  becomes  thicker.  It  becomes  stabilized 
after  60  seconds.  The  calculated  gas  temperatures  are  compared  with 
the  measurements  in  Figure  35.  The  reasons  for  discrepancy  in 
temperature  measurements  at  the  beginning  of  the  test  are  not  clear 
to  us,  We  are  going  to  look  into  this  problem.  Otherwise,  the  cal¬ 
culations  agree  reasonably  well  with  the  measurements.  Gas  tempera¬ 
ture  calculations  for  test  runs  5A  and  14A  are  compared  with  actual 
averaged  temperature  measurements  in  Figure  36  and  37.  The  model 
needs  fine  tuning  to  get  a  better  agreement. 

The  gas  concentrations  of  CO,  HF  and  HCN  are  compared  with 
actual  measurements  in  Figures  38  and  39.  The  calculations  show 
reasonable  agreement  with  the  test  results  ac  the  early  stage  of 
testing.  The  disagreements  become  obvious  after  180  seconds.  This 
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FLAME  SPREAD  RA1ES  AND  DIRECTIONS 


Cylindrical  Flame  volume;  radius 
height,  hf. 
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Thermal  Discontinuity  Position 


Figure  36 


Gas  Temperature  Test  M 


CO  Concentration  Test  1AA 
a — &Avg.  cone.  60  in  above  floor  (±200  ppm) 
-  DACFIR3  upper  zone 
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Figure  38 


HF  and  HCN  Concentration  Test  14A 


Figure  39 


is  possibly  due  to  the  input  material  data  which  gives  constant  species 
release  rate  during  the  test  period, 

The  current  version  of  DACFIR-3  needs  further  refinements  and 
fine  tuning.  The  model  is  physically  sound  and  the  numerical  pro¬ 
cedures  arc  proven  workable  and  economical.  The  major  shortcomings 
of  the  model  validation  are  the  reliability  and  availability  of 
material  properties  as  input  data.  In  particular,  the  autoignition 
data  and  the  flame  spread  rate  data  were  obtained  in  a  laboratory 
scale  apparatus  and  may  not  be  directly  applicable  to  a  real  full- 
scale  fire.  There  is  a  need  to  correlate  the  laboratory  data  to  a 
full-scale  test.  The  species  release  rates,  which  were  obtained  from 
a  laboratory-scale  apparatus,  require  further  examination.  The  effects 
of  reduced  oxygen  concentration  on  spread  and  emission  rates  need  to 
be  incorporated  into  the  model  once  the  data  becomes  available 
(Figure  40) . 

The  additional  refinements  are  shown  in  Figure  41.  The  computer 
code  needs  improvements  and  rearrangement  to  streamline  its  computa¬ 
tions.  In  order  to  account  for  the  radiation  on  vertical  and  ceiling 
surfaces,  the  circular  cylindrical  flame  model  may  not  be  adequate. 

A  final  draft  report  will  be  completed  and  forwarded  to  the  FAA 
for  review  in  two  months.  The  computer  code  and  the  listing  are 
available  through  the  FAA, 
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POTENTIAL  REFINEMENTS  AND  PROBLEM  AREAS 


Figure  40 


ADDITIONAL  REFINEMENTS 


Figure  41 


CORRELATION  WORK  AND  FLAME  SPREAD 


JAMES  QUINTIERE 

Head  of  the  Math  Modeling  Group,  Center 
for  Fire  Research,  National  Bureau  of 
Standards  (NBS) .  Ph.D.  Mechanical 
Engineering,  New  York  University. 
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CORRELATION  WORK  AND  FLAME  SPREAD 


James  Quintiere 

Center  for  Fire  Research 

National  Bureau  of  Standards 

We  have  two  projects  for  which  we  are  responsible  to  the  FAA, 
as  shown  in  Figure  1.  One  is  entitled,  "Correlation  between  Laboratory- 
scale/subscale/full-scale  Fire  Tests."  We  are  a  half  year  into  that 

project.  The  second  is  a  project  on  the  development  of  some  new  con¬ 
cepts  in  flame  spread  methods.  I  will  try  to  outline  what  we  are  up 
to  in  these  two  projects  in  the  following  slides. 

The  output  of  the  correlation  work  will  be  presented  to  the  FAA 
at  the  end  of  the  year.  We  are  in  the  midst  of  it  right  now,  and 
both  myself  and  Bill  Parker  are  involved  in  this.  We  are  focusing 
on  three  elements  in  looking  at  the  relationship  between  test  methods, 
scale  modeling  and  full-scale  fire  results.  Those  elements  might  be 
composed  of  flammability,  burning  rate,  flame  spread,  smoke,  and 

toxicity  (Figures  2  and  3).  The  correlation  work  consists  of  litera-  j 

ture  reviews  in  two  major  areas,  shown  in  Figure  4,  We  want  to  find 

out  what  has  been  done;  specifically,  how  do  test  methods  correlate  \ 

with  full-scale  results,  what  analyses  has  there  been  of  fire  test 

methods  in  the  past,  and  the  same  goes  with  regard  to  scale  modeling.  I 

J 

We  are  excluding  pressure  modeling  in  the  scale  modeling  review.  We 

are  just  looking  at  atmospheric  modeling  techniques  and  how  well  they  ! 

have  performed,  We  are  approaching  this  beyond  detail  and  routine 

features  of  a  literature  review.  We  want  to  see  if  we  can  understand 

the  underlying  features  of  some  of  these  test  methods.  ' 

The  analyses  (Figure  S)  may  call  for  some  generic  mathematical 
modeling  in  simple  terms  of  what  the  test  method  is  trying  to  do. 

We  need  to  get  at  what  the  significant  outputs  of  these  test  methods 

are.  In  this  process,  we  might  be  able  to  identify  what  are  the  more 

important  things  that  are  being  measured  and  have  relevance  as  compared  < 
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Figure 


Figure 


Figure 


to  maybe  what  is  just  empirical  to  provide  a  ranking  order  of  materials 
in  this  testing  apparatus. 

The  objective  (Figure  6)  is  to  relate  this  to  the  FAA's  fire 
scenario  that  they  are  studying;  i.e.,  the  postcrash  fire.  This  is 
the  focus  of  the  output  of  this  review  and  basically  that  is  the  ob¬ 
jective.  Through  this  literature  review,  through  some  analyses  of 
test  methods  and  an  understanding  of  what  the  FAA  is  up  to  in  their 
full-scale  postcrash  fire  tests,  we  hope  to  develop  a  strategy  for 
making  recommendations  on  what  kinds  of  test  methods,  what  kinds  of 
data,  what  kind  of  approach  should  go  into  unravelling  this  and  come 
up  with  a  risk  assessment  for  this  particular  scenario.  This  is  the 
objective.  We  are  in  the  midst  of  this  work  which  will  be  reported 
at  the  end  of  the  project  in  September  or  October  of  this  year. 

The  second  project,  the  flame  spread  test  method  development,  is 
about  to  commence.  We  have  already  done  developmental  work  for  ma¬ 
terials  in  a  room  fire  which  will  be  used  as  guidelines  to  the  approach 
to  this  project.  We  are  preparing  to  start  testing  some  concepts. 

Now,  we  are  putting  together  an  apparatus  to  get  this  project  underway. 

What  is  this  project  all  about?  We  are  attempting  to  develop  two 
concepts  that  will  allow  us  to  predict  mathematical  relationships  for 
rate  of  flame  spread  in  terms  of  measured  quantities  from  small-scale 
test  apparatuses.  We  view  flame  spread  in  a  very  simple  two-element 
mode.  One  is  so-called  creeping  spread,  which  is  spread  against  the 
flow  of  gases,  against  the  wind  if  you  prefer,  shown  in  Figure  7. 

This  is  like  spreading  downward  on  a  wall  or  spreading  laterally  on 
a  wall.  The  other  mode  of  flame  spread  is  wind-aided  flcme  spread, 
also  shown  in  Figure  7.  This  could  be  flame  spreading  up  a  wall  or 
spreading  under  a  ceiling  and  the  wind  can  be  generated  by  the  fire 
itself.  In  this  sense,  we  are  trying  to  separate  the  two  extreme 
modes  of  flame  spread  and  develop  some  test  method  strategy  for  this. 
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Figure 


In  order  to  do  this,  we  would  like  to  explore  materials  that 
are  distinctly  different  to  cover  all  aspects  of  fire  properties  and 
flame  spread.  The  list  of  materials  is  shown  in  Figure  8.  We  would 
choose  at  NBS  three  materials  that  tend  to  represent  what  people  look 
at  in  the  building  side  of  fire  spread.  The  FAA  would  select  three 
materials  that  are  more  relevant  to  the  aircraft  fire  problem.  In 
this  way,  we  would  come  up  with  a  wide  range  of  materials.  Tenta¬ 
tively,  we  have  selected  wood  which  may  be  a  particle  board  or  a 
fiberboard,  and  PMMA  which  is  a  favorite  specimen  for  a  lot  of  people. 

We  would  like  to  produce  some  data  consistent  with  those  from  former 
studies.  Low  density  polyurdthane  foam  has  the  unique  property  of 
being  very  low  density  and  has  interesting  flame  spread  characteristics. 
Panel  material  from  an  aircraft  is  a  very  complex  multilayered  material. 
Seat  cushion  and  perhaps  a  carpet  will  also  be  the  candidate  materials. 
This  is  a  tentative  set  of  materials  for  flame  spread  studies. 

The  approach  to  the  creeping  spread  problem  will  be  outlined  here 
(Figure  9).  We  have  two  test  apparatuses.  One  would  be  operated  to 
study  flame  spread  downward.  The  other  will  be  operated  to  study  flame 
spread  laterally.  Both  are  radiant  panel  type  apparatuses.  A  dis¬ 
tributed  amount  of  radiation  shines  on  the  sample  such  that  the  high 
radiation  flux  is  at  the  end  of  the  sample  ignited;  the  low  radiation 
end  of  it  is  the  direction  toward  which  the  flame  is  spreading.  By 
appropriate  operation  and  analyses  of  the  data,  we  hope  to  derive  a 
relationship  that  would  yield  flame  spread  as  a  function  of  some 
material  properties.  The  apparatus  that  we  have  been  operating  at  the 
present  time  is  in  the  lateral  mode.  The  radiant  panel  is  inclined. 

In  that  orientation,  it  shines  radiation  of  about  5  watts  per  square 
centimeter  at  the  igniting  end  to  about  a  couple  of  tenths  of  a  watt 
down  at  the  far  end.  The  sample  is  about  8"  tall  and  about  2-1/2'  in 
length  ii.  this  lateral  direction.  The  flame  spread  can  be  seen  moving 
or.  thi..  sample.  This  apparatus  was  designed  by  Alex  Robertson  of  NBo. 

It  has  been  tested  up  to  now  for  possible  use  by  the  International 
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Standards  Organization.  They  are  considering  such  an  apparatus,  but 
they  are  not  using  it  in  the  sane  way  that  we  are  planning  for  this 
work.  Work  with  this  apparatus  has  also  been  supported  under  Coast 
Guard  sponsorship.  We  are  preparing  to  initiate  work  here. 

Questions  of  whether  we  need  to  extrapolate  in  some  fashion  to  a 
turbulent  flame  are  unclear  to  us  right  now.  1  think  you  will  agree 
that  if  we  are  considering  flame  spread  down  at  the  leading  edge,  the 
flame  is  going  to  look  the  same  whether  it  is  6'  tall  or  6"  tall. 

Flame  spread  in  the  lateral  direction  may  be  another  question.  We 
can  do  such  things  as  treat  the  boundary  layer  and  make  it  turbulent 
and  look  for  differences  in  the  apparatus.  Testing  a  larger  sample 
with  this  apparatus  is  not  too  practical  at  the  present  time.  This 
is  a  convenient  way  of  getting  a  relationship  by  testing  one  material 
at  one  time.  It  will  yield  flame  velocity  as  i  function  of  flux  or 
more  important  as  a  function  of  surface  temperatures.  This  is  what 
we  are  trying  to  achieve. 

What  we  are  seeking  is,  by  using  this  apparatus  and  by  using 
some  specific  way  of  operating  it  and  interpreting  the  data,  a  re¬ 
lationship  shown  in  Figure  10.  The  results  of  testing  the  material 
will  be  this  parameter  Cf  and  T^  ,  so-called  ignition  temperature  for 
this  mode  of  flame  spread.  We  have  studied  this  and  a  paper  on  this 
subject  that  will  be  coming  out  in  Fire  and  Materials.  Some  flame 
spread  data  are  shown  in  Figures  11  and  12. 

The  approach  for  deriving  wind-aided  flame  spread  rate  is  out¬ 
lined  as  follows.  Flame  spread  upward  or  under  a  ceiling  is  very 
rapid.  Current  techniques  that  are  used  to  judge  the  flammability 
of  materials  in  that  mode  are  not  scientifically  based.  The  challenge 
then  is  how  can  we  make  some  measurements  for  materials  to  obtain  an 
expression  for  this  rapid  spread  upward  or  under  ceiling.  We  doi  't 
believe  that  we  can  achieve  that  by  making  a  measurement  where  we 
watch  the  flame  moving.  We  can  achieve  that  by  making  measurements 
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EXTERNAL  RADIANT  FLUX  (IV,  cm2  J 

Figure  11 


Theory  Material  Exp  Data  Source 
1.  Plywootf-vinyl  O  Beyris  (9) 


SURFACE  TEMPERATURE  [=Cj 

Figure  12 


for  a  fixed  amount  of  material  burning  in  a  so-called  steady  fashion. 
The  flame  is  not  spreading,  but  it  is  burning  and  its  burning  rate 
may  be  changing  with  time.  We  need  to  expose  it  to  radiation.  The 
heat  transfer  characteristics  of  that  flame  above  the  burning  part 
of  the  sample  is  important,  shown  in  Figure  13.  Eventually  in  the 
course  of  this  work,  we  expect  to  develop  an  apparatus  in  which  we 
have  radiation  shining  on  the  material.  At  an  inert  place  above  the 
material,  we  will  measure  flame  height  and  flame  heat  flux.  We  are 
not  quite  sure  how  to  put  all  this  together  in  a  convenient  test 
method  package  yet.  This  is  the  goal  of  this  work. 

In  the  meantime,  before  we  build  a  test  apparatus,  we  would  look 
at  measurement  parameters  (Figure  14)  which  are  effective  parameters 
for  real  materials--heat  of  vaporization.-  heat  of  reaction,  effective 
air- -fuel  ratio,  and  maybe  flame  length  a»\l  heat  transfer.  Bill  Parker 
is  working  on  some  techniques  to  measure  at  least  the  first  three 
quantities.  We  will  measure  them  in  an  apparatus  which  is  known  as 
the  NBS  Rate  of  Heat  Release  Colormeter.  It  has  a  number  of  radiant 
panels  and  can  be  operated  with  a  sample  vertically  or  horizontally. 

We  will  look  at  the  sample  vertically.  We  will  operate  it  in  a  mode 
in  which  we  are  using  the  oxygen  consumption  technique  to  measure 
energy  release  rate.  The  sample  will  be  on  a  load  cell  so  we  will 
measure  the  weight  loss  continuously.  We  will  measure  the  energy 
release  rate  by  oxygen  consumption.  From  that,  we  hope  to  be  able 
to  deduce  these  properties. 

The  analyses  on  the  test  data  are  shown  in  Figure  15.  What  we 
seek  is  to  look  at  the  effect  of  heat  flux.  We  need  to  ccuple  into 
any  flame  spread  results  the  effect  of  time,  Obviously,  a  thick 
material  will  burn  a  longer  time  than  a  thin  material  and  these 
differences  have  to  be  accounted  for  ultimately.  Hopefully,  we  can 
develop  a  flame  spread  lelationship  that  will  functionally  be  written 
down  as  opposed  to  just  symbolically  written  down.  We  don’t  feel 
that  we  at  NBS  have  the  ability  to  generate  all  of  this  work  and  we 
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Figure  14 


are  getting  some  special  analytical  support  that  will  assist  us  to 
develop  a  relationship.  It  may  not  be  a  unique  situation,  but  it 
will  be  a  step  in  the  direction  that  the  modelers  need.  On  the 
other  side  of  this,  it  will  still  be  a  way  in  which  people  can  rank 
mr terials--even  if  they  don't  want  to  use  the  results  of  this  equa¬ 
tion.  Hopefully,  we  will  have  done  it  with  a  better  scientific  basis 
than  people  had  the  resources  to  do  20  to  30  years  ago  when  some  of 
the  current  flame  spread  test  methods,  that  are  currently  in  exis¬ 
tence,  were  developed. 

QUESTION: 

What  is  the  factor? 

JAMES  QUINTIERE: 

The  Cf  factor  has  things  like  thermal  conductivity,  and  heat  transfer 
from  the  flame  in  it.  What  would  be  interesting  is,  if  we  develop 
some  techniques  for  the  upward  flame  spread  and  have  some  techniques 
for  the  downward  flame  spread,  to  see  if  some  of  these  parameters  are 
consistent  between  the  two  techniques.  For  example,  will  ignition 
temperature,  if  derived  from  downward  flame  spread  by  data  analysis, 
be  the  same  as  ignition  temperature  for  upward  flame  spread  that  we 
fit  in  the  model  like  this?  The  same  goes  for  these  other  things, 
the  constants  like  thermal  conductivity, 

QUESTION: 

When  oxygen  consumption  technique  is  used,  do  you  consider  reactions 
as  stoichiometric? 

JAMES  QUINTIERE: 

The  only  thing  you  can  say  about  oxygen  consumption  is  that  you  could 
find  a  lot  of  examples  where  it  looks  like  it  was  a  sound  technique. 
There  may  be  some  that  chemists  can  turn  up  that  don't  work  so  well. 

It  seems  that  from  what  is  in  the  literature  that  you  can't  say  the 
technique  is  going  to  work,  but  works  for  enough  of  the  cases  that 
it  looks  like  it  is  OK. 

QUESTION: 

Is  preheat  level  included  in  the  test  matrix? 

JAMES  QUINTIERE: 

Yes,  the  work  on  the  radiant  panel  test  for  the  lateral  spread-pre¬ 
heating  is  an  important  consideration  in  assuring  that  we  made  the 
proper  analyses  from  tY*  results.  The  reason  is  that  the  rate  of  flame 
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spread  is  not  a  unique  function  of  heat  flux.  It  is  only  a  unique 
function  of  the  heat  flux  if  the  sample  that  you  are  heating  has 
equilibrated  as  a  result  of  that  external  radiant  heat  flux.  That 
time  for  equilibration  is  the  preheat  time  that  we  need.  It  is 
different  for  each  material.  It  is  something  that  we  have  to  fix; 
otherwise  to  use  such  a  technique  as  a  test  method,  the  operator 
of  the  test  would  always  have  to  know  what  that  preheat  time  is  for 
the  material. 

QUESTION: 

How  do  you  measure  heat  of  vaporization? 

JAMES  QUINTIERE: 

To  measure  the  heat  of  vaporization,  it  probably  would  be  best  to  do 
it  in  some  inert  atmosphere,  It  is  not  practical  though  with  the 
apparatus  that  we  are  considering  to  use  right  now.  We  don't  know 
if  we  measure  heat  of  vaporization  with  char  oxidation  whether  that 
is  the  same  e^e  you  would  measure  in  an  inert  atmosphere  or  what  the 
Jifferen:‘\s  ore. 
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UNDSAFE  CODE  APPLIED  TO  AIRCRAFT 
CABIN  FIRE  MODELING 

K.T.  Yang 

University  of  Notre  Dame 

We  have  heard  quite  a  bit  about  the  relationship  between  room 
fire  and  aircraft  cabin  fire.  It  should  be  quite  clear  that  despite 
the  differences  in  the  scenarios  and  also  material  characterizations, 
there  may  still  be  basic  fire  modeling  techniques  applicable  to  both 
situations.  Our  project  at  Notre  Dame  is  also  part  of  the  FAA  math 
modeling  effort  through  the  Interagency  Agreement  between  FAA  and 
NBS.  The  principal  investigators  and  their  associates  are  listed  in 
Figure  1. 

The  objective  of  our  project,  shown  in  Figure  2,  is  to  use  a 
two-dimensional  field  model  (UNDSAFE)  that  we  have  developed  in  the 
last  several  years  and  apply  it  to  an  aircraft  cabin  fire  problem. 

The  specific  things  we  would  like  to  look  at  are  effects  of  fire 
source  strength  and  location.  There  are  several  different  places  in 
a  fuselage  where  a  fire  could  be  initiated.  We  would  also  like  to 
take  a  look  at  the  effect  of  doorway  configuration.  UNDSAFE  is  a 
two-dimensional  model,  The  only  change  we  can  make  is  the  height  of 
a  doorway  opening.  We  would  also  like  to  take  a  look  at  the  effects 
of  seating,  if  seats  would  actually  burn.  Finally,  we  would  like  to 
take  a  look  at  the  effect  of  vertical  venting.  We  have  done  some 
preliminary  work  in  this  particular  area.  It  is  a  very  effective 
way  of  venting  the  combustion  products  out  of  a  room.  We  would  like 
to  take  a  look  at  that  for  aircraft  cabin  fire  venting  problems. 

UNDSAFE  code  was  developed  for  room  fires  and  we  have  since 
made  some  modifications  on  the  basic  code  to  simulate  aircraft  cabin 
fire.  Major  modifications  and  current  progress  are  shown  in  Figure 
3.  The  heat  losses  along  the  ceiling  to  the  outside  vent  become 
important  factors.  On  the  basis  of  some  very  crude  modeling,  we 
can  also  take  into  account  the  additional  heat  release  given  off  by 


127 


COMPUTER  MODELING  OF  AIRCRAFT  CABIN  FIRE  PHENOMENA 


GRANT  NB81NADA  2000 
CENTER  FOR  FIRE  RESEARCH 
NATIONAL  BUREAU  OF  STANDARDS 

TO 

DEPARTMENT  of  AEROSPACE  and  MECHANICAL  ENGINEERING 
UNIVERSITY  of  NOTRE  DAME 


INVESTIGATORS: 

DR.  JOHN  R.  LLOYD 
DR.  A.  MURTY  KANURY 
DR.  K.  T,  YANG 


RESEARCH  ASSOCIATE 
K,  SATOH 


GRADUATE  ASSISTANTS 

X.  Y.  ZHONG 
H.  S.  KOU 


Figure  1 


128 


**■  H*.-  ’r  .  w.  , .  -  . 


OBJECTIVES:  TWO-DIMENSIONAL  FIELD  MODEL  (UNDSAFE) 

SIMULATION  OF  LONGITUDINAL  SPREAD  OF  HOT 
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EFFECTS  OF  FIRE  SOURCE  STRENGTH  AND  LOCATION 
EFFECTS  OF  DOORWAY  CONFIGURATION 
EFFECTS  OF  SEATING  AND  BURNING  SEATS 
EFFECT. OF  VERTICAL  VENTING 


Figure  2 
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Figure  3 


the  seat  when  seat  surface  temperature  reaches  a  pre-set  level. 
Finally,  we  include  an  additional  equation  for  smoke  concentration, 
assuming  the  heat  source  is  also  a  smoke  source.  Smoke  will  be 
propagated  throughout  the  cabin.  Currently  we  are  working  on  two 
separate  problems.  One  is  a  simulation  of  a  full-scale  cabin  fire 
experiment.  The  second  one  is  a  simulation  of  fire  in  a  cabin  with 
seats. 

A  decision  was  made  last  September  in  Dayton  to  use  both  a  zone 
model  by  C.  MacArthur  and  a  field  model  by  Notre  Dame  to  simulate  a 
full-scale  fire  experiment  at  NASA/Johnson  Space  Center.  Test  3B, 
which  was  a  fire  inside  a  737  fuselage  with  seats  and  two  openings 
was  chosen  to  be  modeled  (Figure  4) .  We  are  going  to  make  a  com¬ 
parison  at  the  60-second  point  into  the  fire.  During  the  test,  data 
indicated  that  fuel  weight  loss  rate  was  almost  constant,  as  shown 
in  Figure  5.  This  simplifies  the  situation,  even  though  the  actual 
code  can  actually  incorporate  that  into  the  computation. 

The  second  one  takes  into  consideration  the  heat  losses  through 
the  ceiling.  There  is  a  heat  transfer  coefficient  for  the  fuselage. 
Obviously,  it  takes  some  time  before  the  heat  loss  effect  becomes 
important,  At  the  60-second  point,  we  did  not  feel  that  the  problem 
was  so  serious  that  you  had  to  include  heat  loss  through  a  ceiling. 
The  dimensions  of  a  737  test  article  and  instrument  locations  are 
shown  in  Figures  6,  7,  8,  and  9. 

Figure  6  shows  the  fuselage  configurations  of  a  737.  It  is 
almost  symmetrical  and  looks  like  a  two-dimensional  configuration, 
other  than  the  fact  that  two  doors  are  in  the  aft.  In  order  to  use 
a  two-dimensional  code,  we  have  to  make  some  modification  to  accomo¬ 
date  that. 

When  you  talk  about  a  simulation  of  this  type,  you  really  have 
to  stop  and  think  about  what  you  are  doing.  Because  of  the  many 
parameters  in  this  model  and  also  because  we  use  a  two-dimensional 
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Figure  4 
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TEST  3B  P9N  WEIGHT 
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737  CABIN  MOCK-UP  DIMENSIONS  —  TEST  3B 


Figure 


(/Sc.lc*) 


p*-**  tot 


.  0  -  'T‘hr~+c*u.p/,  i  L  2 C  /W6  /j'u«~*  /  ^-e 

fb  1  c*n  -5  Z+fftJsJ  J 

^1.  .0^  ~  S *<#/»«  /  K-/  l<r  !•*)(*. 

^  CSv*  'rf'ftj 

0*.  CC>  4  *'  *v  J 

{  a.J  J-  u 

!e(<  -  '■'^4/ 

®  “  Cohtt  iVw  liv  £  tfryj 

\J 

C*  f  5*v  »>S  ,  £»<«  ■r<r 

Cvc/jJ  ici«( 

i'-a-'cLtfr  5 

/  ^  “/^Vn/i^  *J  /*  S'^»MJ|'/  Ca/s/,'*. 

✓  /  ' 

^-*'*’3  4— »  arc  .  .< 

•  V‘  <kf'  V£«;t  £«■*•*,  ^  «•*  *»  &~*t  ^ 

Figure  7 

AN  LOCATION  AND  DETAILS  -  TEST  3B 


code,  only  the  two-dimensional  equivalents  of  throe-dimensional 
phenomena  are  simulated.  We  hope  that  we  will  be  able  to  do  this 
because  the  basic  configuration  is  very  close  to  two-dimensional, 
but  there  are  places  in  the  geometry  where  three-dimensional  effects 
become  quite  important.  We  vary  these  parameters  to  get  a  reasonable 
agreement  with  the  experimental  data.  The  basic  equivalents  are  heat 
load,  fire  shape,  and  doorway  heights,  shown  in  Figure  10.  We  do 
have  a  loose  constant  in  a  turbulence  model  which  would  enable  us  to 
employ  different  mixing  levels  to  see  how  that  would  affect  the  re¬ 
sult. 

We  do  not  anticipate  that  a  perfect  agreement  between  simulation 
and  experimental  data  will  be  obtained.  Besides,  there  were  also 
uncertainties  in  experimental  measurements,  as  indicated  in  Figure  11. 
The  numerical  values  of  two-dimensional  equivalent  quantities  for 
heat  loss,  door  height  and  fire  shape  are  listed  in  Figure  12.  The 
two-dimensional  equivalent  heat  load  feels  hotter  (349  KW)  than  the 
actual  experimental  value  (235  KW) .  The  door  height  is  1.05  meters 
compared  to  1.56  meters.  This  is  understandable  because  in  a  three- 
dimensional  case,  an  additional  chocking  effect  occurring  at  the  door¬ 
way  cannot  be  modeled  by  a  two-dimensional  code.  We  have  a  4  x  4 
cell  of  a  fire  source  at  the  bottom  and  2x8  cells  on  the  top  to 
generate  heat.  This  arrangement  will  give  a  ratio  of  height  and  base 
of  a  fire  to  obtain  a  desirable  fire  shape. 

Figure  13  shows  a  comparison  of  calculations  with  experimental 
data.  The  top  portion  of  the  simulation  is  quite  good  throughout 
the  length  of  the  fuselage. 

Figure  14  gives  the  appearance  that  width  is  very  large  compared 
to  fuselage  length.  This  is  actually  not  the  case,  We  plotted  it 
this  way  simply  because  this  is  the  way  that  data  were  obtained.  The 
width  is  small  compared  to  the  length  of  a  fuselage.  Figure  15  shows 
calculated  temperatures  at  four  different  heights  for  a  heat  source 
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RATIONALE  FOR  THE  SIMULATION  STUDY 

TWO-DIMENSIONAL  SIMULATION  MODEL 
2D  vs ,  3D 
DOORWAY  HEIGHT 

DETERMINATION  OF  2-D  EQUIVALENT  OF 

HEAT  LOAD 
FIRE  SHAPE 
DOORWAY  HEIGHT 
CONSTANT  IN  TURBULENCE  MODEL 


Figure  10 


FULL-SCALE  EXPERIMENT  SIMULATION'  EXPECTATIONS 


UNCERTAINTIES  IN  EXPERIMENTAL  MEASUREMENTS 

TEMPERATURE  MEASUREMENTS  CLOSE  TO  FIRE 
EXTENT  OF  FIRE  PLUME 

NO  ATTEMPT  TO  OBTAIN  PERFECT  AGREEMENT 


Figure  11 
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of  349  KW.  The  effects  of  heat  input  rates  (400  KW  and  600  KW)  on 
temperature  profiles  are  shown  in  Figure  16.  The  effect  of  energy 
distribution  is  shown  in  Figure  17. 

The  next  ten  figures  have  to  do  with  the  second  exercise  that 
we  have  gone  through.  We  are  going  to  simulate  fires  in  a  wide-body 
cabin  with  seats.  The  geometrical  arrangement  of  seating  is  modeled 
by  a  two-dimensional  equivalent.  In  this  model,  six  seats  are  set 
along  the  same  line  with  a  heat  source  taking  place  between  the  third 
and  fourth  seats  from  the  left.  Figure  18  shows  the  temperature  pro¬ 
files  in  a  cabin  without  seats.  The  two-zone  effect  is  clearly  demon¬ 
strated. 

Two  distinct  seat  configurations  were  used  in  the  model.  The 
first  seat  configuration  has  a  solid  seat  bottom,  and  the  second  seat 
configuration  has  an  opening  under  a  seat  cushion. 

The  total  input  for  this  particular  computation  was  700  KW  and 
also  when  temperatures  exceed  about  1000°F,  each  cell  will  generate 
an  additional  5  KW.  A  sequence  of  fire  spread  from  an  early  fire  at 
0.96  second  to  a  fully  developed  fire  at  32.15  seconds  is  shown  in 
Figures  19  to  27. 

A  fire  was  first  confined  in  between  the  seats  (Figures  19  and 
20).  Hot  gases  rising  from  the  fire  reach  the  ceiling  and  start  to 
move  along  the  ceiling  (Figures  21  to  25).  At  5.76  seconds  into  the 
fire,  the  hot  gases  reach  the  two  openings  at  both  ends.  Due  to 
different  soffit  heights,  the  flow  patterns  are  different  (Figure  26). 
A  two  layer  effect  is  clearly  indicated.  At  32.15  seconds,  the  fire 
becomes  fully  developed.  The  neighboring  seats  are  heated  and  the 
hot  gases  at  the  top  become  thicker  and  descend  down  to  the  lower 
layer  (Figure  27).  UNDSAFE  code  also  calculated  velocity  vectors 
and  species  concentrations  inside  an  aircraft  cabin.  The  gas  re¬ 
circulations  near  the  openings  are  clearly  demonstrated  by  the 
changing  of  vector  directions  at  the  corners.  This  effect  has  not 
been  simulated  by  a  zone  model  calculation, 
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Figure  25 


We  have  made  a  different  set  of  calculations  with  seats  having 
openings  underneath.  The  calculations  show  that  some  hot  gases  rising 
from  a  fire  are  redirected  to  other  seats  because  the  cool  air  cir¬ 
culates  through  openings  under  the  seats.  We  have  additional  burning 
of  seats  and  additional  heat  is  generated  in  that  particular  area. 

This  is  an  interesting  comparison  between  two  sets  of  calculations. 

It  implies  that  fire  spread  can  be  limited  if  air  circulation  is 
limited.  This  favors  a  blocking  or  partition  concept  in  an  aircraft 
design.  From  a  safety  viewpoint  this  is  not  very  conclusive.  It 
implies  that  a  higher  seat  back  would  have  a  beneficial  effect  as  far 
as  fire  safety  is  concerned.  Another  interesting  point  is  that 
with  an  open  bottom,  temperatures  are  much  lower  than  at  surround¬ 
ing  areas  away  from  the  heat  source  because  circulating  air  cools  the 
flame  temperature.  There  is  a  trade-off.  Further  studies  are  required 
to  clear  this  interesting  problem. 

QUESTION : 

Where  is  the  fire  located? 

K.T.  YANG: 

The  fire  is  located  at  the  center  of  the  fuselage. 

QUESTION: 

The  temperatures  are  low  compared  to  that  in  a  fire,  less  than  200°C? 
K.T.  YANG: 

Yes,  the  temperatures  are  low.  We  were  concerned  about  this.  That 
is  the  reason  why  a  flame  is  shaping  up  like  this  in  the  model. 

QUESTION: 

Everything  shown  in  Figure  14  is  calculated? 

K.T.  YANG: 

Yes,  everything  is  calculation.  If  you  connect  all  the  points  with 
a  straight  line  you  will  get  a  curve.  If  you  connect  the  experimental 
points  by  straight  lines,  you  get  a  different  curve,  which  indicates 
what  actually  is  going  in  the  main  fuselage. 

QUESTION: 

Those  calculations  say  that  the  plume  is  about  20  or  30  centimeters 
wide? 
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K.T.  YANG: 
Yes. 


QUESTION: 

Is  this  because  you  are  forcing  it  into  two  dimensions? 

K.T.  YANG: 

No,  I  don't  think  this  is  the  case.  Don't  forget,  this  scale  is  mis¬ 
leading.  The  fuselage  is  something  like  56  feet. 

QUESTION: 

All  right,  maybe  it  is  40  centimeters  in  width.  For  a  two-dimensional 
model  maybe  that  is  some  justification.  But  for  a  radial  model,  it  is 
clearly  going  to  be  maybe  a  foot  wide  or  more. 

K.T.  YANG: 

Additional  data  will  be  needed  to  determine  that. 

QUESTION: 

The  maximum  temperature  is  only  250°C  over  the  plume  where  combustion 
is  located.  There  is  no  combustion? 

K.T.  YANG: 

No,  the  combustion  occurs  in  the  plane. 

QUESTION: 

How  do  you  define  the  plane?  The  lowest  level  of  the  temperature 
should  be  much  higher. 

K.T.  YANG: 

I  think  what  you  are  getting  at  is  some  skepticism  on  the  part  of 
people  who  have  run  fire  tests  and  made  measurements  inside  planes, 
and  the  skeptical  position  is  with  regard  to  the  possibility  that 
you  might  measure  the  temperatures  which  are  not  higher  than  250°C 
in  a  place  where  there  is  burning  going  on. 

We  have  often  thought  about  how  really  accurate  measurements  are  and 
we  can  get  the  temperatures  much  lower  than  they  are  if  you  take  into 
account  the  radiation  factors.  You  have  to  be  very  careful  about 
exactly  the  measurement  conditions. 

The  problem  might  be  that  there  was  a  very  coarse  thermocouple 
grid  in  the  experiment.  Instead  of  having  the  thermocouples  on  the 
axis  of  a  fire,  it  may  be  one  foot  off.  If  you  are  trying  to  match 
the  numerical  values  of  two  temperatures,  then  you  are  way  off. 

This  is  t ho  host  data  we  have.  We  have  to  have  some  way  to  make 
some  comparison  just  to  see  what  kind  of  equipment  and  data  we  are 
t  a  1 k i ng  about . 
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MODELING  HEAT  FLUXES  FOR  AIRCRAFT 


Ronald  Alpert 
Factory  Mutual  Research 

and  Engineering  Corporation 

The  title  of  this  project  is  computer  modeling  of  aircraft  cabin 
fire  phenomena. 

We  are  going  to  formulate  a  few  efficient  computer  subroutines 
that  could  be  used  in  a  comprehensive  zone  model.  I  am  going  to 
describe  the  plans  for  this  project. 

The  first  task,  shown  in  Figure  1,  is  to  develop  integral  models 
of  fire  spread  under  corridor  ceilings.  The  integral  models  can  be 
very  efficient  on  computer  time  and  yet  reasonably  accurate.  The 
geometry  in  Figure  1  is  this  one  where  a  flow  exists  along  the  wall 
and  the  ceiling.  The  wall  will  be  combustible,  but  the  ceiling  may 
or  may  not  be  combustible.  The  side  walls  are  to  confine  the  flow 
at  the  wall  and  the  ceiling.  The  plane  view  on  top  shows  what  might 
happen  if  the  ceiling  is  combustible. 

A  flame  occurs  and  the  flame  front  progresses  down  the  ceiling. 
That  is  the  general  view  of  what  we  are  looking  at. 

Factory  Mutual  is  under  an  FAA  contract  to  conduct  an  experimental 
study  on  physical  modeling.  We  have  run  intermediate-scale  and  small- 
scale  experiments  on  ceiling  burning.  A  good  deal  of  data  exists 
from  these  experiments  which  could  be  used  as  a  comparison  with 
theoretical  prediction  calculated  from  integral  models. 

Figure  2  delineates  the  specific  objectives  of  Task  1  work. 

During  the  first  year,  we  will  be  looking  at  the  first  two  topics. 
First,  we  will  want  to  validate  existing  integral  models  of  combus 
tion  in  fire  plumes.  Second,  we  will  also  develop  and  validate 
integral  models  for  wall  fires. 

There  are  several  different  types  of  integral  models  that  we 
have  developed  at  Factory  Mutual  for  fire  plume  combustion.  The 
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TASK  1: INTEGRAL  MODELS  OF  FIRE 
SPREAD  UNDER  CORRIDOR  CEILINGS 


Validation  of  Integral  Combustion 
Models  for  Turbulent  Fire  Plumes 
and  Wall  Fires 

Formulation  of  an  Integral  Model 
for  Reacting,  Turbulent  Wall- 
Ceiling  Flows 

Solution  of  the  Ceiling  Flow 
Combustion  Model  with  Comparison 
to  FMRC-FAA  Experiments 

Formulation  of  Transient,  Under- 
Ceiling  Fire  Spread  for 
Incorporation  into  Zone  Models 

Figure  2 
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first  model  was  developed  by  Dr.  Francesco  Tamanini.  His  integral 
model  is  a  modification  of  the  numerical  techniques  for  reaction  in 
a  buoyant  turbulent  plume.  He  makes  assumptions  in  order  to  simplify 
his  model  and  develop  a  rather  efficient  integral  model  for  buoyant 
turbulent  combustion  in  a  fire  plume.  This  is  one  model  that  is 
quite  promising  for  use  with  fire  plumes.  This  integral  model  will 
allow  low  cost  predictions  on  the  rate  of  burning  in  fire  plumes. 

We  need  experimental  data  for  verification  of  this  model. 

Another  model  of  buoyant  combustion  in  the  plume  is  Dr.  John 
de  Ris'  stochastic  mixing  model,  involving  the  evolution  of  probability 
density  function  in  the  plume.  Again,  we  have  a  simplified  model  which 
requires  comparison.  Optimization  from  experiments  in  this  process  is 
actively  being  pursued  right  now. 

An  experimental  apparatus  developed  by  Dr.  F.  Tamanini  was  used 
to  verify  the  integral  models  of  plume  combustion.  It  has  a  water- 
cooled  chamber  with  gas  burners.  By  raising  the  burner  up,  various 
levels  of  a  plume  can  be  experimentally  studied.  The  flame  enters  a 
duct  with  gas  analyses  instruments  downstream.  We  can  determine  the 
chemical  composition  of  the  products  and  degree  of  reaction  at  that 
level  in  the  plume.  In  addition,  a  radiometer  is  mounted  in  the  side 
of  the  chamber  so  we  can  look  it  a  th  .n  slice  of  flame.  The  radiant 
output  from  the  slice  of  flame  can  '  compared  with  the  heat  release 
at  that  same  level.  A  typical  result  is  that  the  energy  release 
rate  integrated  in  the  plume  is  a  function  of  height  above  the  burner. 
It  can  determine  the  fraction  of  the  fuel  converted  to  carbon  monoxide. 
For  the  first  time,  we  have  some  hard  data  on  where  in  the  fire  plume 
the  chemical  reaction  is  occurring.  This  same  apparatus  will  be  used 
by  Mike  Delichatsios  to  obtain  measurements  of  air  entrained  in  the 
plume--a  technique  very  similar  to  that  used  by  Professor  Ed  Zukoski 
from  California  Institute  of  Technology.  This  one  apparatus  will 
allow  us  to  make  these  critical  comparisons  between  theory  and  experi¬ 
ment,  to  validate  models  of  plume  combustion.  Once  this  validation 
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has  taken  place,  we  can  go  on  to  looking  at  the  problem  of  the  wall 
fire  and  developing  integral  models  for  wall  fires. 

The  remainder  of  the  first  year  we  will  be  looking  at  formula¬ 
ting  integral  models  for  the  wall/ceiling  combustion  configuration 
and  extending  the  wall  fire  integral  model  to  combustion  under  the 
ceiling.  In  the  second  year  of  the  program,  we  hope  to  solve  the 
ceiling  flow  combustion  model  and  compare  predictions  from  the  theory 
with  the  experimental  data.  Finally,  once  we  have  the  steady  solution, 
we  will  formulate  a  transient  under  ceiling  fire  spread  model  by  con¬ 
sidering  the  transient  case  to  be  just  a  succession  of  steady  burning 
situations.  The  integral  solution  could  then  be  incorporated  in 
existing  zone  models. 

Task  2  of  this  NBS  grant  deals  with  the  three-dimensional  solu¬ 
tion  of  fire  heat  transfer  in  an  aircraft  cabin.  The  situations  we 
will  be  looking  at  are  the  radiant  heat  transfer  from  a  pool  fire 
outside  the  aircraft  to  the  interior  of  the  cabin  where  some  penetra¬ 
tion  occurs  (Figure  3).  This  is  under  quiescent  wind  conditions. 

With  the  outside  pool  fire  and  entrained  air  from  the  cabin,  the  flame 
has  been  drawn  into  the  upper  part  of  the  cabin.  The  flame  penetrates 
down  the  cabin  and  forms  a  hot  ceiling  layer  going  down  the  length 
of  the  cabin.  We  will  be  looking  at  the  situation  where  we  have 
flame  penetration  into  a  cabin,  looking  at  the  radiant  flux  and  con¬ 
vective  flux  to  arbitrary  targets  within  the  cabin. 

Figure  4  shows  the  specific  objectives  of  Task  2.  In  the  first 
year,  we  do  the  first  two  subtasks  and  in  the  second  year  the  last 
two  subtasks.  The  first  subtask  is  to  calculate  radiant  heat  trans¬ 
fer  from  external  pool  fires  to  arbitrary  targets  within  the  air¬ 
craft.  We  will  develop  both  numerical  solutions  and  approximate 
analytical  solutions  so  that  we  can  judge  the  accuracy  of  the 
approximate  solution.  These  solutions  will  be  in  terms  of  parametric 
properties  of  the  outside  pool  fires.  In  the  remainder  of  the  first 
year,  we  will  be  estimating  heat  transfer  due  to  flame  penetration. 
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TASK  2:  THREE  DIMENSIONAL  SOLUTION  FOR 
FIRE  HEAT  TRANSFER  IN  AN  AIRCRAFT  CABIN 

Radiant  Heat  Transfer  from  an 
External  Pool  Fire  to  an  Arbitrary 
Target  inside  an  Aircraft 

Estimation  of  Heat  Transfer 
due  to  Flame  Penetration 

Improved  Calculation  of 
Penetrating  Flame  Heat  Transfer 
with  Results  from  TASK  1 

Computer  Subprogram  for  Efficient 
Calculation  of  Heat  Transfer 
Rates  from  Reacting  Wall~Ceiling 
or  Plume~Cei 1 ing  Flows  to  an 
Arbitrary  Target  within  the  Aircraft 

Figure  4 
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The  assumed  geometries  or  thicknesses  of  the  penetrating  flame  and 
the  properties  of  the  penetrating  flame  will  be  used  to  calculate 
the  convective  and  radiant  heat  transfers  to  the  ceiling  and  arbitrary 
targets  in  the  aircraft.  In  the  second  year,  based  on  calculations 
made  in  Task  1,  we  will  be  looking  at  ceiling  layer  combustion.  We 
will  predict  properties  of  the  layer,  i.e.,  thicknesses  of  the  layer 
and  temperatures.  We  will  use  that  information  for  improving  the 
calculation  of  heat  transfer  to  targets  within  the  aircraft. 

Finally,  in  the  remainder  of  the  second  year,  we  will  try  to 
develop  efficient  subroutines  for  the  calculation  of  heat  transfer 
from  either  the  penetrating  flame  or  the  ceiling  layer  flame  or 
combustion  products  in  the  aircraft  cabin  to  arbitrary  products 
within  the  aircraft.  This  is  our  plan  for  this  project. 

We  are  just  beginning  and  Mike  Delichatsios  has  started  work 
on  the  wall  fire  combustion  problem  and  has  made  some  real  progress 
there. 

QUESTION: 

Are  you  going  to  attempt  to  work  out  a  method  of  inserting  the 
results  that  you  get  into  a  zone  model? 

RONALD  ALPERT: 

It  would  be  very  nice  if  we  could  do  that.  It  depends  on  timing. 

If  we  have  something  developed  on  time  during  the  contract  period, 

I  think  we  would  look  at  it.  We  have  the  capability,  for  instance, 
for  running  the  Harvard  program  at  Factory  Mutual. 

QUESTION: 

Why  did  you  pick  this  particular  geometry  for  the  hulli’  You  left 
out  all  the  return  flow  problems.  I  realize  that  it  makes  life 
easier,  but  in  a  real  hull,  it  might  not  be  totally  unimportant. 

There  is  some  data  that  has  been  obtained  that  shows  that  you  can 
actually  get  local  air  built  up  near  the  fire  which  would  cause 
a  couple  of  lengths  of  the  hull  to  be  completely  afoul  with  smoke. 
There  are  problems  of  that  sort  that  are  associated  with  returning 
flow.  Are  you  only  interested  in  the  very  thin  layer  on  the  top? 
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RONALD  ALPERT: 


We  wanted  to  tackle  a  problem  and  solve  to  a  degree  that  we  really 
believe  the  answer.  We  don't  want  to  go  out  further  than  we  think 
we  can  catch.  We  had  enough  problems  with  getting  a  combustion  model 
working  correctly.  We  had  enough  challenges  with  this  one  without 
taking  on  further  challenges  at  this  time. 

QUESTION: 

It  seems  to  me  you  would  want  to  address  the  problem  of  a  pool  fire 
outside  the  doorway.  You  are  working  on  a  piece  of  that.  I  wonder 
if  you  are  going  to  address  or  consider,  even  when  there  is  no  wind 
acting  on  the  fire  plume,  there  is  some  sporadic  intermittent  pene¬ 
tration  of  that  plume  into  the  cabins.  Have  you  made  any  considera¬ 
tion  or  will  you  address  in  your  work  or  will  you  hope  others  will 
address  some  work  on  how  to  describe  that  phenomena- -how  do  you  see 
that  related  to  your  problem. 

RONALD  ALPERT: 

I  hope  someone  will  describe  that  phenomena.  I  don't  see  that  we 
are  going  to  predict  a  random  penetration  of  the  plume.  We  may  simu¬ 
late  it  by  saying  we  have  a  wall  fire  and  a  respectively  black  body 
source  there  or  a  wall  fire  on  one  side  and  then  have  a  ceiling 
flow  generated  by  that  wall  fire  or  have  some  assumed  type  of  plume 
being  a  fraction  of  the  pool  fire. 
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ENCLOSURE  MODELS  APPLIED  TO  AIRCRAFT 


Henri  Mitler 
Harvard  University 

My  discussion  is  going  to  very  briefly  explain  the  current  status 
of  the  Harvard  Computer  Fire  Code  and  how  this  is  applicable  to  the 
ceiling  jet  problems  of  the  FAA  effort  in  general. 

The  Harvard  Computer  Fire  Code  is  a  deterministic  model  of 
fires  burning  in  enclosures.  At  the  moment,  the  enclosures  that 
we  are  considering  are  rectangular  (i.e.,  a  room)  which  has  a  number 
of  vents  (limited  to  five).  We  can  also  handle  the  behaviors  of  up 
to  five  objects,  at  least  one  of  which  is  burning  and  the  rest  are 
to  be  considered  targets.  The  floor  will  be  considered  an  object. 

The  math  model  is  deterministic.  The  computer  program  is  modular 
so  that  we  can  remove  a  subroutine  if  we  wish  and  substitute  an  even 
simpler  one  or  a  better  one  or  a  more  complicated  one.  This  includes 
not  only  physical  subroutines  but  also  numerical  subroutines.  We 
have  basically  two  numerical  subroutines.  One  is  a  successive  sub¬ 
stitution  method  of  solving  an  enormous  set  of  simultaneous  equations. 
Another  is,  as  C,  MacArthur  pointed  out,  a  Newton-Raphson  technique. 

The  fires  we  are  looking  at  are  pre-flashover  fires.  Harvard 
Computer  Fire  Code  Version  Five  of  this  model  is  about  to  come  off 
the  drawing  board.  A  tape  of  version  five  should  be  made  soon  and 
if  anyone  is  interested  in  obtaining  that  tape,  please  see  me.  This 
model  can  handle  several  types  of  fires.  One  is  a  growing  fire,  such 
as  igniting  a  piece  of  polyurethane  foam  or  a  mattress  or  anything 
else  and  watching  this  thing  grow.  Moreover,  a  fire  can  be  set 
initially  and  can  be  ignited  at  some  point  down  the  road  either  by 
autoignition  due  to  the  charring  surface  of  the  flammable  target 
reaching  an  ignition  temperature  or  by  contact  with  the  flame.  We 
model  flames,  vent  flow  rates,  plume,  species  production,  the  spread 
rate  for  a  growing  fire,  convective  heat  transfer,  etc. 
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Figure  1  shows  a  schematic  of  the  enclosure  fire.  A  flame  is 
modeled  by  a  cone  of  hot  gas  which  is  assumed  to  be  a  grey  emitter 
with  uniform  properties.  The  flame  temperature  is  chosen  by  the 
user.  We  modeled  it  with  an  optical  absorption  coefficient  of  1.55 
reciprocal  meter.  We  also  have  a  plume  model  for  the  hot  gases  rising 
from  the  flame.  We  used  the  Morton-Taylor-Turner  point  source  plume 
model  shown  in  Figure  2.  A  virtual  point  source  is  located  below  the 
fire  surface  and  the  plume  itself  will  assume  either  a  top  hat  model, 
which  is  actually  what  I  use  now,  or  a  Gaussian  profile.  It  makes 
very  little  difference.  The  radius  of  the  plume  is  in  effect  the 
radius  of  the  fire  at  the  burning  base.  A  virtual  part  of  the  plume 
is  below  the  fire  base  and  the  real  part  of  the  plume  has  an  air  en¬ 
trainment  coefficient  which  is  assumed  constant.  Nevertheless,  in 
spite  of  the  simplicity,  the  plume  model  has  worked  quite  adequately. 

The  flow  rate  of  hot  gas  and  air  are  shown  in  Figure  1.  The 
flame  and  plume  go  up  to  the  ceiling  and  form  a  hot  layer  which  gets 
deep  in  time;  then  buoyancy  carries  it  out.  We  solved  the  ventilation 
equation  effectively  the  same  way  as  C.  MacArthur  pointed  out  to  you. 

In  fact,  we  have  drawn  independently  exactly  the  same  basic  equations. 

We  model  the  species  concentration  by  assuming  that  any  particular 
object  gives  rise  to  a  constant  mass  fraction  of  carbon  monoxide, 
carbon  dioxide  and  smoke  which  consists  of  mostly  soot  and  possibly 
hydrocarbon.  One  source  of  these  mass  fraction  data  is  burning  poly¬ 
urethane  foam  by  Tewarson.  I  used  the  numbers  that  he  developed.  It 
is  a  very  simple  approximation  with  single  numbers.  Nevertheless,  the 
results  are  reasonable  for  all  the  species  except  carbon  monoxide. 

We  have  to  use  experimental  results  for  the  flame  spread  rates. 

We  could  not  get  the  spread  rate  from  first  principles.  It  is  possible 
to  have  an  expression  which  gets  the  correct  spread  rate  for  the  open 
flame,  but  then  corrects  for  the  effect  of  the  feedback  radiation  from 
the  hot  layer,  hot  ceiling,  the  walls,  etc.  Again,  we  get  quite  rea¬ 
sonable  results  by  doing  that. 
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We  also  deal  with  all  the  radiation  exchanges  between  surfaces 
of  hot  gases  and  plume.  The  radiative  transfer  calculations  are 
fairly  reasonable.  Convective  heating  and  convective  cooling  are 
calculated  by  using  a  heat  transfer  coefficient. 

The  standard  ^oom,  which  is  a  room  used  by  Factory  Mutual  Re¬ 
search  Corporation,  is  8'  x  12'  x  8'.  The  burning  rate  is  a  function 
of  time  and  is  very  well  predicted.  The  radius  of  the  fire  is  a  func¬ 
tion  of  time.  The  calculated  and  experimental  data  are  compared  and 
shown  in  Figure  3.  The  surface  temperature  stays  almost  constant 
until  approximately  2-1/2  to  3  minutes  into  the  burn.  The  calcula¬ 
tion  results  and  the  experimental  data  are  shown  in  Figures  4  and  5. 
The  temperature  of  the  gas  layer  reaches  a  peak  of  about  900°F  and 
then  goes  off.  I  have  not  yet  worked  with  a  model  that  has  any  kind 
of  extinguisher. 

QUESTION: 

What  is  the  oxygen  concentration? 

H.  MITLER: 

The  minumum  that  I  find  is  on  the  order  of  6  percent.  It  is  bounded 
by  a  fraction  of  4  from  ambient  oxygen  and  for  carbon  dioxide;  the 
experiment  does  something  like  this--it  goes  up  expeditiously,  throws 
some  heat  and  falls  off. 

The  oxygen  fraction  measured  data  is  compared  with  calculated 
results  and  both  curves  are  shown  in  Figure  6.  The  calculated  value 
starts  at  the  ambient  value  and  follows  the  same  pattern  as  the 
experimental  data. 

The  calculated  and  experimental  results  on  CO,  CO^  and  smoke 
production  are  shown  in  Figure  7.  The  calculated  CO  doe;  not  match 
the  measurement.  The  measured  CO  starts  out  at  a  much  lower  value 
and  increases  to  a  factor  of  5  or  more  higher  than  what  is  calculated. 
It  may  well  be  that  this  is  the  different  nature  of  the  burning  of 
vitiated  air  or  non-vitiated  air.  Very  likely  it  is  because  of  the 
distribution  of  flaming  combustion  products,  particularly  carbon 
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monoxide  in  the  flame  itself.  At  the  moment,  I  have  not  made  any 
attempts  to  "correct”  the  model  for  that  because  it  is  a  less 
serious  thing.  The  CO  at  the  moment  is  the  worst  predictor  we  have 
made  so  far. 

Absorption  calculations  are  shown  in  Figure  8.  Two  curves  of 
infrared  absorption  coefficient  for  the  layer  are  obtained  by  using 
two  different  subroutines.  One  is  extremely  cimple  by  assuming  that 
the  absorption  coefficient  is  proportional  to  the  soot  concentration. 
The  curve  ABSRB-2  starts  low,  reaches  a  peak  and  then  falls  off. 

The  other  curve,  ABSRB-3,  which  is  practically  on  top  of  ABSRB-2, 
comes  from  a  subroutine  which  is  rather  more  elaborate  and  takes  into 
account  the  carbon  dioxide,  1^0 ,  and  soot  concentrations.  Both 
curves  can  be  fitted  equally  well  to  experimental  data.  It  also 
indicates,  but  is  not  obvious  from  these  curves,  that  you  cannot  use 
the  soot  concentration  in  the  layer  to  get  this  result  unless  you 
specifically  take  into  account  the  soot  deposition  on  the  ceiling 
and  walls.  Fifteen  to  20  percent  of  the  soot  must  be  lost  by  deposi¬ 
tion.  The  absorption  coefficient  would  be  too  large  by  that  amount 
all  the  way  through. 

I  have  estimated  the  computer  CPU  time  required  for  a  run  of 
1,000  seconds  of  fire  time.  The  CPU  time  for  the  version  four  of 
this  computer  code  which  would  be  1  -  1-1/2  years  old  is  220  seconds 
on  a  PDP-VAX  machine,  a  large  mini-computer.  This  was  done  about 
eight  months  ago,  April  1980.  Thirty  percent  computer  time  could  be 
saved  by  changing  the  numerical  algorithms  and  the  logical  flow  of 
the  numerics  but  maintaining  the  same  physics.  What  I  did  here  was  to 
make  a  substantive  substitution  algorithm  more  robust.  At  least  for 
the  standard  run,  it  never  has  invoked  the  Newton-Raphson  method. 

One  iteration  of  this  algorithm  takes  20  milliseconds,  but  one 
iteration  of  the  Newton  algorithm  takes  600  milliseconds.  The 
present  version  has  improved  and  expanded  physics.  Last  summer  and 
fall  we  improved  the  logical  flow  of  the  program  and  made  it  completely 
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sound.  The  trimmed  computer  program  economizes  the  calculations  and 
CPU  time. 

The  computer  program  is  quite  flexible  and  has  b'.en  designed 
for  the  convenience  of  the  user.  We  have  run  various  kinds  of  cases. 
We  have  calculated  fire  behaviors  in  a  very  long  corridor  quite  suc¬ 
cessfully.  We  have  run  a  case  where  gasoline  was  spilled  on  a  carpet. 
In  the  computer  program,  the  gasoline  vapor  was  treated  as  one  object, 
the  saturated  carpet  was  treated  as  a  second  object,  and  the  unsatu¬ 
rated  carpet  as  a  third  object.  Again,  we  had  quite  reasonable  re¬ 
sults.  We  have  run  a  variety  of  other  carpet  fire  cases,  and  the 
numerical  results  occasionally  ran  into  a  snag.  However,  in  version 
six  we  have  had  to  try  several  new  algorithms,  one  of  which  is  by 
Dr.  de  Ris.  We  eventually  hope  to  have  a  set  of  numerical  programs 
which  are  sufficiently  powerful  that  practically  any  problem  could 
be  solved. 

We  intend  in  this  version  six  to  expand  major  actions.  We  do 
not  yet  have  wail  fires,  the  effect  of  fire  extinguishers,  and  multi¬ 
ple  fires. 

As  far  as  the  applications  to  the  FAA  efforts,  we  are  going  to 
expand  the  model  capabilities.  Professor  Emmons  has  already  written 
down  some  equations  to  model  the  ceiling  jet  and  this  will  most 
likely  be  incorporated  into  the  model. 

QUESTION : 

What  are  the  details  of  your  new  successive  substitution?  The  old 
one  seemed  to  give  you  a  lot  of  problems. 

H.  MITLER: 

There  are  essentially  two  changes  I  made.  The  original  successive 
substitution  mode  was  a  Jacobian  technique  where  we  used  the  old 
value  to  calculate  all  the  new  values.  The  first  change,  we  went  to 
a  Gauss-Seidel  method  which  means  that  any  time  we  calculate  the 
variable,  we  use  that  updated  value  instead  of  the  previous  iteration 
value.  The  second  part  which  made  it  much  more  robust  is  that  after 
re-alternating  Gauss-Seidel  iterations,  I  then  took  the  mean  value 
of  the  last  iteration  and  the  current  iteration  and  that  seemed  to 


make  it  much  more  robust  and  it  works  very  well.  Another  part  of 
the  numerical  package  is  that  I  try  to  avoid  using  the  Newton  tech¬ 
nique  as  much  as  I  can  because  of  the  large  amount  of  work  involved 
in  solving  the  Jacobian.  I  use  the  Newton  technique  when  necessary 
and  then  I  immediately  shift  out  of  that  to  a  grid-size  where  I 
don’t  change  the  Jacobian. 

QUESTION: 

Have  you  used  double  precision  on  the  VAX? 

H.  MITLER: 

Yes  and  no.  I  have  not.  John  Randall,  a  graduate  student  working 
with  us,  has  used  double  precision.  He  finds  that  sometimes  double 
precision  is  needed  in  solving  the  Jacobians  and  has  to  go  with  double 
precision. 

QUESTION: 

You  make  some  comparisons  between  the  computed  upper  layer  tempera¬ 
ture  and  the  experimental  upper  layer  temperatures  that  should  be 
interfaced.  How  do  you  compute  this  from  your  experimental  data? 

H.  MITLER: 

There  are  three  racks  of  thermocouples  (front,  middle  and  rear). 

Each  rack  has  a  dozen  or  more  thermocouples.  We  weighed  the  numbers 
in  some  reasonable  way  for  all  those  temperatures  and  took  an  average 
of  those  numbers.  We  are  also  working  on  a  couple  of  equations  which 
tie  them  together. 
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THERMOCHEMICAL  MODELING  OF 
BURNING  AIRCRAFT  MATERIALS 

Kumar  Ramohalli 

Jet  Propulsion  Laboratory 

This  is  quite  a  different  kind  of  work  from  what  you  have  heard 
until  now.  I  have  decided  to  spend  a  little  time  explaining  what  we 
have  done  so  far  in  proper  perspective. 

A  plan  for  this  presentation  is  outlined  in  Figure  1.  I  will 
show  the  basic  approach  and  objectives  that  we  are  taking.  I  would 
also  like  to  give  some  background  information  on  what  we  have  done 
so  far,  since  many  of  you  have  not  seen  our  work  before.  Next,  I  shall 
talk  about  our  work  with  the  FAA  Technical  Center  during  FY'80. 

Finally,  our  future  plan  for  FY'81  and  beyond  on  thermochemical 
modeling . 

The  objectives  (Figure  2)  of  thermochemical  modeling  are  as 
follows: 

•  Predict  fire  and  smoke  behavior  using  only  the 
ingredient  thermochemical  property  values  and 
the  geometry  employed,  preferably  in  a  non- 
empirical  fashion. 

•  Suggest  economical  methods  for  better  materials 
and  transfer  these  methods  to  industry. 

•  Progressive  steps  in  complexity. 

Figure  3  shows  the  relationship  between  thermochemical  modeling, 
fire  modeling,  and  material  research.  The  fire  modeling  covers  all 
reactions,  heat  and  mass  transfers  in  the  gas  phase.  The  interface 
between  fire  modeling  and  thermochemical  modeling  occurs  at  the 
boundary  layer  of  the  material  surface  where  char  reradiation, 
oxidative  degradation,  porosity  and  heterogenous  effects  are  accounted 
for.  For  the  condensed  phase,  subsurface  degradation,  polymer  frag¬ 
mentation,  radicals  diffusion,  charring,  outgasing,  and  layered 
structural  members  are  modeled.  It  is  assumed  that  one  gram  of  con¬ 
densed  material  becomes  one  gram  of  vapor  phase  material  without 
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any  drastic  change  during  the  phase  transition.  For  the  material  re¬ 
search,  we  model  chemical  reaction  in  simple  terms  of  Arrhenius  factor 
and  activation  energy.  The  hydrogenation  kinetics,  bond  strength, 
agglomeration  and  nucleation  kinetics  are  also  considered.  The  actual 
species  concentrations  at  the  material  surface  play  an  important  role 
in  the  modeling. 

The  oxidative  degradation  process  and  fragment  size  vaporizing 
are  shown  in  Figure  4.  The  undisturbed  polymer  has  very  large  molecu¬ 
lar  weights,  which  is  sometimes  in  the  range  of  several  million.  As 
the  distance  to  the  surface  decreases,  the  temperature  is  increased 
and  the  polymer  starts  to  break  down.  It  finally  leaves  the  surface 
but  the  vaporizing  surface-fragment  size  is  not  clearly  known.  There 
is  rarely  any  analytical  treatment  on  surface-fragment  size  vaporizing 
(FSV)  and  the  assumptions  of  monomers,  dimers,  and  trimers  have  yet 
to  be  experimentally  verified.  It  is  modeled  that  FSV  is  related  to 
surface  temperature  and  vapor  pressure  and  a  Clausius-Clapeyron  type 
of  vapor  pressure  equilibrium  is  used.  Data  on  vapor  pressure  from 
the  American  Petroleum  Institute  on  hydrocarbons  are  shown  in  Figure 
5  and  Figure  6.  It  is  shown  that  the  vapor  pressure  is  a  function  of 
temperature  for  various  hydrocarbon  species  and  there  is  a  definite 
trend  that  vapor  pressure  increases  with  surface  temperature  for  a 
given  hydrocarbon  molecular  weight. 

The  equations  for  subsurface  degradation  are  shown  in  Figure  7. 
The  first  equation  is  a  one-dimensional  energy  balance  equation  near 
the  surface.  The  second  equation  is  the  polymer  degradation  equation. 
The  boundary  conditions  are  specified  for  temperatures  and  molecular 
weights  at  the  surface  where  a  fire  exists  and  at  a  location  far  away 
from  the  fire,  The  equations  are  highly  nonlinear  but  an  analytical 
solution  is  possible  because  the  exponential  factor  is  greater  than 
30.  Some  details  of  solution  for  charring  materials  are  shown  in 
Figure  8.  The  energy  equation  is  defined  for  both  char  and  solid 
regions  with  matched  heat  flux  boundary  conditions.  A  second  order 
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Figure 


ordinary  differential  equation  is  non-dimensionalized,  Simultaneous 
solutions  for  burning  rate  of  material  and  fragment  size  at  char  layer 
are  obtained  by  iterative  techniques.  The  temperature  at  the  con¬ 
densed  phase  is  determined  from  heat  flux.  The  effect  of  temperature 
of  the  material  on  burning  rate  is  clearly  indicated  by  the  solution. 
The  solution  also  indicates  that  there  is  a  temperature  gradient  in¬ 
side  the  char  and  the  material  undergoes  a  drastic  change  in  chemical 
composition. 

The  specific  tasks  in  FY'80  for  the  FAA  are  listed  in  Figure  9. 
Carpet  and  aircraft  seat  cushion  are  the  test  materials.  The  objec¬ 
tives  are  to  predict  burning  behavior  of  these  materials  under  the 
following  conditions: 

1.  Incident  radiation 

2.  Self-sustained  flame 

3.  Treatment  with  flame  retardant 

4.  Thermal/physical  thickness 

5.  Spalding  "B"  number 

6.  Ambient  pressure 

In  this  presentation,  I  have  time  for  only  one  of  the  tasks. 

The  burning  rate,  or  the  reaction  rate,  of  a  seat  cushion  with 
flame  retardant  is  shown  in  Figure  10.  A  thermal  gravimetric  analysis 
(TGA)  technique  was  used  to  measure  the  temperatures.  In  a  burning 
case,  the  actual  hydrocarbon  concentration  at  the  surface  can  be  quite 
small  but  still  measurable.  It  is  in  the  order  of  one  or  two  percent . 
The  reaction  rate  is  higher  in  a  nitrogen-filled  environment  than  that 
in  an  air-filled  environment.  Figure  11  shows  the  reaction  rates  of 
carpet  which  is  used  on  an  aircraft  floor  deduced  from  TGA  technique. 

The  predicted  burning  rates  of  carpet  are  shown  in  Figure  12.  It 
is  predicted  that  the  heat  flux  supplied  by  the  flame  to  the  surface 
is  not  sufficient  to  sustain  the  flame  at  high  burning  rates.  Heat 
flux  from  other  sources  are  required  to  supplement  the  radiative  flux 
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PREDICTED  BURNING  RATES  OF  CARPET 


Figure  12 
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from  the  material  flame  for  sustaining  the  combustion.  The  material 
surface  under  radiative  flux  could  raise  the  surface  temperature  and 
become  self-ignited.  The  burning  rate  for  self-ignition  is  low. 

The  flame  actually  starts  moving  farther  and  farther  out  as  the 
burning  rate  goes  up  and  heat  flux  from  the  flame  decreases.  If  no 
outside  heat  flux  is  supplied,  the  burning  rate  will  be  diminished 
and  extinguished. 

The  specific  tasks  in  FY'81  for  the  FAA  are  summarized  in  Figure 
13.  The  material  to  be  tested  is  multilayered  polymer  such  as  honey¬ 
combed  panels,  polyurethane  foam-neoprene  blocking  layer  with  wool 
and  nylon  as  seat  covers.  We  plan  to  extend  the  analytical  thermo- 
chemical  model  to  a  multilayered  system  to  predict  burning  behavior 
under  various  heat  flux  conditions.  Different  sizes  of  layer  thick¬ 
ness  will  be  tested  to  obtain  an  optimum  combination  of  multiple 
layers.  More  experimental  work  is  also  planned.  A  NBS  smoke  density 
test  chamber  will  be  used  to  compare  the  model  predictions.  The 
material  samples  will  be  tested  under  varying  incident  radiative  flux 
and  the  weight  loss  and  temperature  profile  in  the  sample  material 
will  be  recorded.  The  experimental  data  will  be  compared  with  data 
from  NASA/ ARC.  It  is  expected  that  cooperation  with  other  fire  re¬ 
searchers  will  produce  satisfactory  results. 
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ENCLOSURE  FIRE  DYNAMICS  MODEL 
FOR  INTERIOR  CABIN  FIRES 

Josette  Bellan 

Jet  Propulsion  Laboratory 

First  of  all,  I  would  like  to  say  that  we  are  not  sponsored  by 
the  FAA  or  the  NSF.  Our  sponsor  is  NASA  Headquarters.  I  would  like 
to  talk  about  our  Enclosure  Fire  Dynamics  Model  for  Cabin  Fire  and 
Combustion  Modeling.  The  staff  on  this  project  is  shown  in  Figure  1. 

Why  are  we  studying  this  problem?  Our  motivation  is  listed  in 
Figure  2.  We  want  a  long-term  development  of  a  mathematical  tool 
to  predict  the  progress  or  burning,  temperature,  and  gas  species 
distributions.  Our  objectives  are  to  build  a  model  for  the  pre¬ 
dictability  of  aircraft  fire  characteristics  and  a  dynamic  response 
of  materials  in  an  accidental  fire  environment. 

Tho  plan  of  the  presentation  is  outlined  in  Figure  3.  First  of 
all,  I  am  going  to  give  you  some  background  the  formulation, 
progress,  status  and  finite  elements  numerical  procedures  requirements. 

The  physical  picture  is  shown  in  Figure  4.  It  is  a  section  of 
the  aircraft  cabin  and  fireproof  floors  with  an  entrance  and  an  exit. 
Our  approximations  of  the  practical  situation  are  shown  in  the  lower 
part  of  Figure  4. 

Even  with  these  approximations,  we  have  some  difficulties  in 
modeling  them.  The  difficulties  encountered  in  establishing  a  de¬ 
tailed  fire  model  enclosure  are  listed  in  Figure  5.  We  realize  that 
there  is  an  inefficiency  in  trying  to  approximate  both  wall  and  core 
phenomena  because  of  constraints  of  money,  time,  anu  computer  time 
and  also  the  lack  of  thermophysical  and  thermochemical  constants  for 
various  aircraft  materials.  Our  equations  are  the  unsteady-state 
conservation  equations,  the  field  equations,  which  are  similar  to 
those  used  in  the  Notre  Dame  math  model.  The  general  form  consists 
of  mass,  momentum,  energy,  and  species  equations  shown  in  Figures  6 
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Figure 


difficulties  in  establishing  a  detailed  enclosure  fire  model 
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CONSERVATION  EQUATIONS 
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Figure  6 


and  7.  We  assume  that  Lewis  number  and  Prandtl  number  equal  to  one, 
and  one- step  chemical  reactions  represent  the  complex  fire  chemistry - 
With  these  approximations,  we  can  eliminate  from  all  but  one  equation 
the  reaction  term  so  that  our  equations  arc  simple  to  solve. 

The  boundary  conditions  are  listed  in  Figure  8.  The  inert  wall 
conditions  determine  the  species  and  the  energy  counts  at  the  surface. 
At  the  entrance  we  have  forced  ventilation  and  all  the  species  and 
temperatures  are  prescribed.  At  the  exit  we  are  going  to  compute 
temperatures  and  species  by  forward  extrapolation.  Finally,  we  are 
assuming  in  the  first  stage  of  the  model  that  the  pool  fire  is  burning 
and  wind  velocity  is  zero  in  axial  direction.  The  fuel  evaporates  in 
the  pool.  I  want  to  point  out  that  the  equation  we  are  using  here  has 
a  transient  operation.  It  has  an  evaporation  equation  rather  than 
equilibrium  equation.  It  has  been  found  that  there  are  important 
discrepancies  between  models  that  use  conventional  thermodynamic 
equilibrium  and  this  type  of  equation  for  this  time  duration. 

Finally  we  have  here  a  boundary  equation  that  gives  us  the  energy 
balance  and  the  surface  of  a  pool  fire.  Again,  we  make  a  thin  wall 
approximation  that  relates  to  latency  of  the  evaporation  and  enthalpy 
that  evolve  from  the  surface  in  the  gas  flow. 

The  turbulent  transport  modeling  term  is  shown  in  Figure  9.  The 
equations  that  I  have  given  you  previously  are  correct  equations; 
however,  they  don't  isolate  turbulent  transport.  In  order  to  isolate 
turbulent  transport,  it  is  a  well-known  procedure  that  all  the  de¬ 
pendent  variables  are  expressed  as  a  sum  of  mean  gradients  plus  tur¬ 
bulence.  A  solution  for  the  mean  values  is  sought  and  the  correla¬ 
tion  terms  are  modeled.  These  are  practical  eddy-diffusion  types  of 
models  with  all  the  density  variations  neglected.  The  laminar  trans¬ 
port  variations  are  increased  by  turbulent  contributions.  The  Lewis 
and  Prandtl  numbers  are  assumed  equal  to  one.  It  is  sufficient  to 
specify  only  one  of  the  transport  terms.  W<^  chose  to  do  it  for  the 
diffusivity  because  we  do  have  an  analogy  of  a  turbulent  jet. 
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TURBULENT  TRANSPORT  -  MODELING 


*  WfOTH  Of  THE  POOL 
-  KtlGHT  Of  ENCLOSURE 


how,  we  are  coming  to  a  very  important  part  which  is  the  tur¬ 
bulence  combustion  shown  in  Figure  10.  It  is  a  very  controversial 

subject  and  that  is  why  we  decided  that  in  our  model  we  are  not  going 

to  specify  whether  combustion  rate  is  either  controlled  by  kinetics 
or  by  diffusion  alone.  We  are  going  to  have  to  choose  one  of  two 
processes  depending  on  which  one  is  a  slower  process. 

For  the  kinetic  one  we  have  a  practical  one-step  reaction  model 
and  for  the  diffusion  one  we  have  a  reaction  proportional  to  the 
quantities  that  are  defined  here.  They  are  the  mean  square  of  fuel, 
pure  oxygen  mass  and  enthalpy.  Our  definition  of  diffusivity  and 
length  scale  are  also  defined  for  the  diffusion  controlled  process. 
Going  back  to  our  computation  equations  we  can  write  the  equations 

for  the  g's.  The  problem  can  be  solved  and  in  order  to  solve  it 

we  need  additional  modeling.  In  order  to  find  an  easy  way,  we  are 
making  the  assumption  of  local  equilibrium  of  the  flow  which  means 
that  transient  convection  and  diffusion  tens  are  going  to  be  small 
in  respect  to  the  production  and  dissipation.  We  can  then  solve  the 
equation  in  the  right  form.  I  don't  want  to  go  any  further  than 
that  except  to  point  out  that  contributions  of  different  tens  are 
involved. 

For  the  quantities  that  are  related  to  the  mass  fraction,  we 
have  production  due  to  turbulence  transport  which  is  divided  by  the 
dissipation  due  to  turbulence  and  sink  due  to  combustion.  I  would 
like  to  point  out  that  the  combustion  terms  have  not  been  modeled 
before  and  we  are  going  to  compare  the  calculations  with  data.  The 
enthalpy  equation  has  in  the  numerator  terms  for  turbulence  trans¬ 
port  and  buoyancy,  and  in  the  denominator  terms  for  turbulence, 
combustion,  radiation  and  pressure  effect. 

The  description  of  radiation  model  is  shown  in  Figure  11.  Radia¬ 
tion  in  a  turbulent  flow  is  a  very  important  thing.  In  order  to  model 
radiation,  we  find  the  solution  of  the  intensity  equation  and  assume 
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RADIATION  IN  A  TURBULENT  FLOW 
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attenuation  of  radiation  by  gas  u  small.  It  is  also  assumed  that 

the  emissive  power  of  the  enclosure  suifsce  is  negligible  with  respect 

to  that  of  the  gases,  particularly  the  flame.  A  classical  radiation 
4 

equation  of  T  is  used.  For  radiation  in  a  turbulent  flow,  this 
classical  equation  is  modified  by  an  additional  quantity  which  is 
proportional  to  the  mean  square  temperature.  We  expect  this  quantity 
to  be  very  important.  Radiation  to  surfaces  for  an  opaque  surface 
is  a  practical  relationship.  But  in  order  to  find  the  radiancy  for 
nonisothermal  surfaces,  the  emissive  power  of  the  surfaces  is  neglected. 
We  have  .ays  and  in  each  ray  we  have  segments.  Each  segment  has  a  con¬ 
stant  temperature  and  these  segments  are  summed  wO  give  the  total  radi¬ 
ation  to  the  surfaces.  Therefore,  we  have  a  nonisothermal  radiation 
model . 

The  uncoupling  and  quasi  linearization  of  the  equations  are 
shown  in  Figure  12.  A  general  type  non-linear  operator  is  defined 
for  the  differential  equations.  It  is  very  difficult  to  solve  a  set 
of  non-linear  partial  differential  equations.  A  two-step  approach  is 
used  to  reduce  these  equations  in  a  solvable  form.  The  first  step  is 
to  uncouple  the  dependent  variables.  All  thjr  values  of  dependent 
variables  from  the  previous  time-step  are  considered  known  quantities. 
The  second  step  is  to  linearize  the  equations  around  these  last  time- 
step  quantities. 

A  finite  element  method  that  we  are  using  for  solving  the  equa¬ 
tions  is  outlined  in  Figures  13  and  14.  There  is  a  misconception 
about  the  linear  finite  element  method.  It  is  not  a  finite  difference 

method,  because  it  does  not  need  the  approximation  at  every  single 
point  on  the  grid.  A  finite  element  method  is  not  a  zone  model  which 
deals  with  an  extremely  large  amount  of  spaej.  The  finite  element 
method  which  we  are  using  is  not  a  finite  element  method  which  has 
been  applied  to  the  structural  problems.  However,  there  are  a  number 
of  publications  which  have  been  published  during  the  last  decade  on 
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UNCOUPLING  ANb  QUASILINEARIZATION  Or  THE  EQUATIONS 
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Figure  12 
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the  finite  element  method.  I  would  like  to  explain  to  you  really 
what  it  means. 

A  two  dimensional  cabin  with  "x”  and  "y"  coordinates  is  shown 
in  Figure  13.  A  vertical  axis  is  called  function  "4*".  For  each 
point  on  a  grid,  I  define  a  function  that  looks  like  a  pyramid  with 
two  segments  both  in  the  "x"  and  "y"  directions.  (I  define  a  variable 
in  such  a  manner  that  it  is  the  sum  of  terms  that  is  a  function  of 
time  only  multiplied  by  a  set  of  base  functions  "4h",  "4y  that  are 
functions  of  "x"  and  "y"  only.)  Now,  for  each  equation  I  fou  d  what 
I  called  the  residual  by  using  a  well  documented  Galerkin's  method. 

The  way  I  got  the  residual  is  as  follows: 

The  expression  for  the  coefficients  was  plugged  back  into  the 
uncoupled  quasi  linearized  equations.  What  you  find  is  that  this 
function  is  not  an  exact  measurement  of  the  equation.  What  you  obtain 
is  something  that  I  call  Rr  which  is  an  error  function.  It  is  re¬ 
quired  that  the  error  vector  be  perpendicular  to  the  vector  "<f>"  to 
minimize  the  error.  There  are  more  details  on  finite  elements 
method  in  Figures  14  and  15, 

This  finite  element  method  was  tested  with  a  simple  heat  con¬ 
duction  equation  and  the  results  are  shown  in  Figure  16.  Finally, 
this  method  is  quite  a  bit  more  obvious  than  the  finite  difference 
method.  There  have  been  comparisons  done  on  that  simply  by  solving 
a  convection-diffusion  equation.  The  second  reason  is  that  this 
method  is  very  easily  amendable  to  a  variable  mesh  (variable  in 
space,  variable  in  time  as  well),  which  may  save  computer  space  and 
time. 

The  future  plan  for  the  JPL  fire  math  modeling  effort  is  pre¬ 
sented  in  Figure  17.  The  finite  elements  method  will  be  further 
tested  with  a  non-combustion  flow  problem.  The  program  has  to  be 
coded  and  the  criterion  for  convergence  needs  to  be  determined.  It 
is  necessary  to  have  access  to  a  very  fast  computer  at  NASA  Langley 
Research  Center  to  test  the  program. 
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APPLICATION  OF  THE  FINITE  ELEMENT  METHOD 
TO  A  HEAT  CONDUCTION  EQUATION 


EQ:  4r  '  d(-^4  +  -^4]  ;  D  •  ~  thermal  diffusivity 
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nr  11  +  Jl  it  -  t  t  ■  e  ( h  -  heat  transfer  coefficient 

BV'5  k  a  ;  |  rT  -  normal  to  the  surface  In  the  outward  direction 
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APPENDIX  A 


FAA/NBS  WORKSHOP  ON  MATHEMATICAL  FIRE  MODELING 
March  24-27,  1981 


AGENDA 


March  24,  1981 
8:30  -  8:45 

8:45  -  9:00 

9:00  -  9:30 

9:30  -  12:00 


12:00  -  12:30 
12:30  -  1:30 
1:30  -  1:50 

2:00  -  2:20 

2:30  -  2:50 


3:00  -  3:30 


Technical  Building  Auditorium 
Opening  Remarks 

Wayne  Howell,  Chief  of  Fire  Safety  Branch 
FAA  Technical  Center 

Aircraft  Fire  Scenarios 

Gus  Sarkos,  Fire  Safety  Branch 

FAA  Technical  Center 

FAA  Modeling  Efforts 

Thor  Eklund,  Fire  Safety  Branch 

FAA  Technical  Center 

DACFIR  Model  Workshop 
Charles  MacArthur,  University  of  Dayton 
Research  Institute 
Dayton,  Ohio 

Coffee  Break 

Discussion 

Lunch 

Correlation  Work  and  Flame  Spread 

James  Quintiere,  National  Bureau  of  Standards 

Washington,  D.C. 

UNDSAFE  Applied  to  Aircraft 
K.T.  Yang,  Notre  Dame  University 
South  Bend,  Indiana 

Modeling  Heat  Fluxes  for  Aircraft 
Ronald  Alpert,  Factory  Mutual  Research 
Corporation 
Norwood,  Massachusetts 

Coffee  Break 
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3:30  -  3:50 

Enclosure  Models  Applied  to  Aircraft 

Henri  Mitler,  Harvard  University 

Cambridge,  Massachusetts 

4:00  -  4:20 

Thermochemical  Modeling  of  Burning  Aircraft 
Materials 

Kumar  Ramohalli,  Jet  Propulsion  Laboratory 
Pasadena,  California 

4:30  -  5:00 

Enclosure  Fire  Dynamics  Model  for  Interior 
Cabin  Fires 

Josette  Be) ’an,  Jet  Propulsion  Laboratory 
Pasadena,  California 

5:00  -  5:30 

Discussion 

March  25,  1981 

Director's  Conference  Room 

4th  Floor 

Technical  Building 

9:00  -  12:30 

Tours  of  FAA  Fire  Test  Facilities 

12:30  -  1:30 

Lunch 

1:30  -  5:00 

Ad  Hoc  Committee  (Plumes) 

March  26,  1981 

Director's  Conference  Room 

4th  Floor 

Technical  Building 

9:00  -  5:00 

Ad  Hoc  Committee  (Plumes) 

March  27,  1981 

Director's  Conference  Room 

4th  Floor 

Technical  Building 


9:00  -  5:00  Ad  Hoc  Committee  (Smoke  Movement) 
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